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ABSTRACT 

A new htting methodology is presented which is equally well suited for 
the estimation of low-, medium-, and high-degree mode parameters from m- 
averaged solar oscillation power spectra of widely differing spectral resolution. 
This method, which we call the “Windowed, MuLTiple-Peak, averaged spec¬ 
trum”, or WMLTP Method, constructs a theoretical prohle by convolving the 
weighted sum of the prohles of the modes appearing in the htting box with the 
power spectrum of the window function of the observing run using weights from 
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a leakage matrix that takes into account both observational and physical effects, 
such as the distortion of modes by solar latitudinal differential rotation. We 
demonstrate that the WMLTP Method makes substantial improvements in the 
inferences of the properties of the solar oscillations in comparison with a previous 
method that employed a single profile to represent each spectral peak. We also 
present an inversion for the internal solar structure which is based upon 6,366 
modes that we have computed using the WMLTP method on the 66-day long 
2010 SOHO/MDI Dynamics Run. To improve both the numerical stability and 
reliability of the inversion we developed a new procedure for the identification 
and correction of outliers in a frequency data set. We present evidence for a 
pronounced departure of the sound speed in the outer half of the solar convec¬ 
tion zone and in the subsurface shear layer from the radial sound speed profile 
contained in Model S of Christensen-Dalsgaard and his collaborators that existed 
in the rising phase of Solar Cycle 24 during mid-2010. 

Subject headings: Sun: helioseismology — Sun: oscillations — methods: data 
analysis — methods: numerical 


1. Introduction 


Helioseismology provides a unique opportunity to investigate in great detail the internal 
structure and rotation of the Sun. The starting point for the study of t he solar interior using 
helioseismology c an be identified with the observational confirmation bv iDeubnerl (119751) . and 
independently bv iRhodes et al.l (119771) . of the standing-wave nature of t he solar five-in inute 
oscillations that were observed in th e solar photosphere as proposed by lUlrichl (jl97nl ). and 
independently bv iLeibacher Sz SteinI (Il97lh. The remarka ble qualitative agreement of the 
observations of Deubnei (j 19751) and of Rhodes et ah (1977) with the predictions of the theo¬ 
retical studies meant that the five-minute oscillations could be regarded as a superposition of 
acoustic normal modes that are trapped in the interior of the Sun. However, the frequencies 
of the observed ridges of power in the dispersion plane were syst ematically lower by about 
5% than the theoreti c al pre dictions. This discrepancy allowed iRhodes et ahl (119761) . and 
independently iGoughI (119771) to provide observational estimates of the depth of the solar 
convection zone. These estimates are now believed to be the first helioseismic inferences of 
solar internal structure. 


The observed modes are predominantly p-modes, for which pressure provides the domi¬ 
nant restoring force. Also observed is the /-mode, which at high spheric al ha rmonic degree 
I has the character of a surface gravity wave. Both the IDeubnerl (119751) and IRhodes et ah 
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(1977) studies employed intermediate- and high-degree /- and p-modes, while IClaverie et ah 
fll979l) employed low-degree p-modes that penetrated into the solar core. Today, the line 
of demarcation between low- and intermediate-degree modes is at / = 4 (the highest degree 
that can be observed in integrated light), while that between intermediate- and high-degree 
modes is where individual modes are no longer resolved, or at about / = 300 for the /-mode, 
between I = 200 and / = 150 for the n = 1 through n = A ridges (and at lower values of I for 
the higher radial orders n). 


The observations by iDeubnerl (119751) . [Rhodes et ahl fjl977l) . and IClaverie et ahl (119791) 
raised considerable debate over th e proper characterization of the oscillation modes that 
had been observed. IDeubnerl (jl977l) cited a pre-publication version of iDeubner et ahl (Il979l) 
to claim that the solar p-mode oscillations were truly global in nature, but subsequently, 
Ulrich et ahl (119791) disagreed with this conclusion and argued that “...the modes of greatest 
interest are not globally coherent because of t he eff ect o f convective m otions associated 
with the supergranulation”. Furthermore, iHilll (119801) and iGoughI (Il980li po i nted out that 
the coherence times of up to nine hours that were cited by iDeubner et ahl (119791) and by 
Claverie et ahl (Il980l) were not sufficient to demonstrate the global nature of these modes. 
Today, the term “global helioseismology” refers either to studies that employ the low- and 
intermediate-degree p-modes whose lifetimes are truly long enough for them to be globally- 
coherent or to studies that employ spherical harmonic decompositions that are computed 
from nearly the entire visible solar hemisphere. Studies that do not use modes which have 
such long lifetimes or which are computed from observations that cover much smaller portions 
of the visible hemisphere are considered to employ the tools of “local helioseismology”. 

Depending on the frequency and degree /, the modes propagate within different acoustic 
cavities inside the Sun between two turning points. Outside the acoustic cavity in which 
a mode is propagating, the mode is evanescent. Therefore, the mode characteristics (e.g., 
frequency) are relatively insensitive to conditions outside the associated acoustic cavity, par¬ 
ticularly far from the turning points. Moreover, modes which propagate in the direction of 
solar rotation have higher frequencies than modes with the same resonant properties prop¬ 
agating in the opposite direction. This effect is called “rotational frequency splitting”. The 
amount of splitting depends on the rotation rate inside the acoustic cavity in which the 
mode is propagating as well as on the azimuthal order m of the mode. Utilizing the dif¬ 
ferential penetration and the frequency splitting of the modes allow s the internal structure 
and r otation of the Sun to be inferred, as a function of position (cf. IChristensen-Dalsgaard 
20021) . This possibility of carrying out inversions of the observed frequencies and frequency 
splittings is a key issue in the applications of helioseismology. So far, however, the vast ma¬ 
jority of inversions performed have only included freque ncies and frequency splittings of the 


low- and the intermediate-degree oscillations (see, e.g.. iGough et al.lll996l: [Thompson et ah 
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19961: lAntia fc Chitrelll998l: iKosovichev et al.lll998l: ISchou et al.lll998l) . On the other hand, 
high-degree /- and p-modes have an immense potential in the helioseismic probing of the 
sub-snrface layers of the Snn. This is demonstrated here in Fignredl where the dependence 
of the inner tnrning-point radius on spherical harmonic degree I is shown for three differ¬ 
ent frequencies, spanning the range of the observations. For small I, the inner turning point 
is rather close to the center of the Sun, whereas for higher degrees it moves closer to the 
surface. In particular, for I > 150 the modes are essentially trapped in the outer 65 Mm 
below the solar surface. Therefore, accurate measurements of high-degree mode frequencies 
and splittings allow us to improve our inf erences regarding the large - scale s t ructure and dy¬ 
nami cs of the sub-surface layers (see, e.g., 

2nn2h . 


Rabello-Soares et al.ll2000L l2008al: iDi Mauro et ah 


Due to the use of a modal concept in global helioseismology, the diagnostic potential 
of the data is necessarily limited. Specihcally, to hrst order the standing acoustic modes 
sense only the longitudinally averaged, north-south symmetric average of the internal strat- 
ihcation of the Sun. Moreover, in contrast to solar differential rotation, flows in merid- 
i onal planes (meridional circu l ation ) have only a tiny effect on global oscillation frequencie s 


(IWoodardI I2OOOI: iRoth fc Stixl I2OO8I: ISchad et al.l I2OIII: IVorontsovI I2OIII: I Schad et al.l l2013l) , 


which severely hampers any attempt to detect such flows by global mode frequency analysis. 
Such limitations are avoided in local helioseismology, which is based upon the assumption 
that the solar oscillations locally behave as propagating acoustic wa ves that are scattere d and 
absorbed by local inhomogeneities and advected by local flow helds. iBraun et al.l (119871) were 
the hrst to demonstrate the utility of this approach by showing that propagating acoustic 
waves could be absorbed by the strong magnetic held associated with sunspots, thus po¬ 
tentially providing information about the magnetic held itself. Subsequently, three methods 
of analyzing propagating aco ustic waves in a localized area on th e solar surface were de¬ 


(Duvall et al. 

I 993 I Kosovichev 

1 996bl: 

Gizon & Birch 

2005: Zhao 

2004, 

2008). and acoustic 

holography ( 

)jindsev & Braun 

2 OO 0 I 2 

004). The application of these methods has led to 

spectacular results, as has been demonstrated by, e.g.. 

Giles et al. 

(1997 

). Kosovichev et al. 

(2000). 

Braun & Lindsev (2001 

L), 

Beck et al. J 2 OO 2 ). 

daber et al.l 

(2002 

). For a review of 


local helioseismology we refer the reader to iGizon et al.l (120101) . 


While in recent years the progress made in local helioseismology has been substantial, it 
has been much slower in global helioseismology. This is in part a consequence of the difficul¬ 
ties inherent in the generation of reliable high-degree mode parameters, but misconceptions 
as to the roles of local and global helioseismology may have contributed as well. One such 


misconception is that t he glob a 


mode measurements can be replaced with measurements 


from local analyses. As iReiterl (120071) has demonstrated, at large-scales the local measure- 
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ments are much less precise than the global measurements and, therefore, there is a strong 
complementarity between the local and global techniques. 

The estimation of high-degree mode parameters is made difficult due to the fact that 
high-degree modes cannot be observed as sharp, isolated peaks but only as ridges of power 
comprised of overlapping modes. Because of the asymmetrical distribution of the amplitudes 
of the modes that blend together, the central frequency of each ridge deviates from the 
frequency of the target mode. Hence, to recover the underlying mode frequency from htting 
the ridge, an accurate model of the ridge power as a function of frequency is required. With 
such an accurate model the global analysis provides the most robust estimates of the mean 
structure and rotation of the Sun which are important for testing theories of stellar structure, 
evolution, and differential rotation. 


We began to delve into the htting of solar oscillation spectra in the late-1980s with 
the development of our hrst-generation, or Single-Peak, Averaged-Spectrum, htting method, 
which we will refer to in the following as either Method 1 or the SPAS method, and which 
is briehy described here in Appendix Because of insurmountable problems with the de¬ 
termination of unbiased ridge-ht frequencies we had to abandon this method, however. We 
therefore began to develop our second-generation, or Windowed, MuLTiple-Peak, averaged- 
spectrum, htting method, which we will refer to in the following as either Method 2 or the 
WMLTP method. This method is equally well suited for the unbiased estimation of low-, 
medium-, and high-degree mode parameters from m-averaged solar oscillation power spectra. 
A detailed description of this method will be presented in Section 01 after we have given an 
outline of the data analysis generally employed in global helioseismology in Section [2l and 
after we have addressed the problems inherent in the analysis of high-degree power spectra 
in Section [3l The issue of the sensitivity of Method 2 in terms of both the m-averaging 
procedure and the effective leakage matrix is addressed in Section |5l Sample results from 
Method 2 are presented in Section O In Section [H] we will also demonstrate the substan¬ 
tial improvements that Method 2 makes in the frequencies, linewidths, and amplitudes that 
we generated with this method by comparing them with corresponding quantities that we 
generated using our Method 1. Finally, in Section [7] we hrstly describe a new procedure 
for the identihcation and correction of outliers in frequency data sets that are to be used 
for solar structure inversions, before we present a new structural inversion from a set of 
frequencies computed frorn^^ 66-d av long spectra obtained with the Michelson Doppler Im- 
ager (MDD (IScherrer et al.lll995l) on board the SOlar and Heliospheric Observatory (SOHO) 
(iDomingo et al.lll9951) at the beginning of solar cycle 24 in 2010, followed by our concluding 
remarks in Section [SI 


At this point we note that the current paper is the hrst of a series of three papers. 
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While the focus of the paper at hand is on Method 2, we will present in the second paper 
our third-generation, or Multiple-Peak, Tesseral-Spectrum, fitting method, which we will 
refer to in the following as either Method 3 or the MPTS method. This method directly fits 
the tesseral, zonal, and sectoral spectra at each degree rather than resorting to m-averaged 
spectra as is the case with both Method 1 and 2. In that paper we will also intercompare 
the results obtained from Method 2 and 3, and we will investigate the systematic effects 
introduced in Method 2 by the m-averaging procedure. The purpose of the third paper will 
be the intercomparison of our results obtained from both Method 2 and 3 with results from 
“established” fitting methodologies at low, intermediate, and high degrees. 


2. Data analysis in global helioseismology 

In global helioseismic studies the data reduction typically includes the following major 
steps. First, the observed Dopplergrams V{9,(f),t) are spatially decomposed into spherical 
harmonic coefficients i.e., 

aUt)= [ w{e,cp)v{e,cp,t)Yr{e,cp)da, (i) 

Jv 

where 6 is co-latitude, 0 is longitude, is a spherical harmonic function of degree / 

and azimuthal order m, W{9, 0) is an apodization chosen to reduce the contribution from 
the noise close to the solar limb, "D is the visible hemisphere of the Sun or a portion thereof, 
and t is time. In this step of the data reduction spatial side-lobes are introduced into the 
target spectrum (/, m) because the spherical harmonic functions Y/^{9, 0) are not orthogonal 
on P. However, even if the entire surface of the Sun could be observed spatial side-lobes 
would be present because of the distortion of the mode eigenfunctions by solar latitudinal 
differential rotation, and also because of velocity projection effects. The integral in equation 
([T]) can be very expensive to compute. A typical approach is to use interpolation to remap, 
for given time t, the product W{9,(j))V{9,(j),t) onto some coordinate system in which the 
integration over 0 may be represented as a Fourier transform, so that Fast Fourier Transform 
(FFT) techniques may be applied. In the second step a spectral analysis is carried out for 
each spherical harmonic coefficient ai^rn{t), i.e., 

/ + 00 

ai,m{t)exp{27rwt)dt, (2) 

-oo 

where u = u/27i is cyclic frequency. In this step temporal side-lobes are introduced into the 
target spectrum (/, m) if periodic gaps are present in the time series of the spherical harmonic 
coefficients due to the day-night cycle, say. In practice, the integral in equation ([2]) 
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is first approximated by a discrete Fourier transform over a finite interval of time, which 
then is efficiently computed using a FFT technique. In the third step the power spectrum 
is calculated for each (/, m). The hnal step in the data reduction consists 
in the peak fitting of the power spectra Alternatively, the peaks in the complex 

spectra can be fitted as well. For example, such approach is employed in the htting 

methodology of ISchoul (119921) . 


2.1. Generation of nn-averaged power spectra 


The results presented in this investigation are based upon four different sets of un¬ 
averaged power spectra that were created from observations obtained with the MDI 

instrument during 1996, 2001 and 2010. The MDI was operated on the SOHO spacecraft 
between April 1996 and April 2011. The MDI obse rvations which we h ave employed were 
all obtained as part of the MDI Full-Disk Program fjScherrer et al.lll995l ). and are listed in 
Table [H where we also indicate the naming convention we use in this paper to refer to each 
observing run. 

Time series of Dopplergrams that resulted from each of the four observing runs listed 
in Table [U were converted into / + 1 complex time series (purely real for m = 0) which 
were gap-hlled using an auto-regressive gap Filing procedure based upon the approach of 
Fahlman fc Ulrvchl fll982[) . using a reduction pipeline that was developed at Stanford Uni¬ 
versity for the processing of the MDI data. The resulting gap-filled duty cycles are listed 
in the last column of Table [B Using standard FFT techniques the gap-hlled I + 1 complex 
time series were converted into a group of 2/ -|- 1 zonal, tesseral, and sectoral power spectra 
for 0 < / < 1000 for each of the four observing runs. In doing so, the positive frequency part 
is identihed with m < 0, while the negative frequency part is identihed with m > 0. 


Within each of these four groups of un-averaged power spectra, each target spectrum 
(/, m) contains a number of frequency bins equal to one-half of the number of samples that is 
listed in the corresponding row of column 3 in Tabled] For the MDI instrument, these fre¬ 
quency bins span the frequency range of zero to the temporal Nyquist frequency of 8333 /iHz. 
Within each group of spectra the zonal (i.e., m = 0) spectrum for a given degree, /, contains 
a variable number of sets of isolated peaks (at low- and intermediate degrees) or a set of 
ridges (at higher degrees) of power that correspond to a collection of /- and p-modes. For 
the cases in which the peaks are isolated, each set of peaks consists of a peak for the target 
mode the set of temporal sidelobes, and a set of spatial sidelobes that have leaked 

into the target spectrum from nearby spectra. We will refer to the entire collection of 2/ -|- 1 
target peaks and their spatial and temporal sidelobes that share a common n-value as the 











(n, 1) multiplet. When we refer to the mode {n, 1), we are actually referring to the m-average 
of the 2 / + 1 modes that share the same values of n and 1. 


For the tesseral and sectoral spectra the corresponding peaks in each spectrum are 
shifted to lower frequencies by solar rotation for the spectra having m < 0 , while they are 


shifte d to higher frequencies for the spectra having m > 0 (cf. IChristensen-Dalsgaard et al. 

2nnnh . 


2.2. Procedures for the generation of m-averaged power spectra 

For a given degree, /, the m-averaged power spectrum is dehned as 


(3) 


where the symbol ()m means averaging over m, is computed in a three-step procedure. In the 
hrst step of this procedure, the frequency shift resulting from the effect of solar rotation and 
asphericity is calculated, f or each of the 2/-I-1 un-averaged spectra, using an iterative 

cross-correlation method f Brown 1985 : Tomczvk 1988 : Korzennik 1990) • In this method the 
frequencie s within a multiple t (n, /) are approximated with a polynomial expansion similar 
to that of Duvall et ah f 19861) . i.e.. 


+ F ^ 4”’'^ Pfc (m/L). 


(4) 


k=l 


Here, i>n,i is the frequency of the multiplet (n,/), = l{l + 1), are the so-called 

frequency-splitting coefficients, and Pk is the Legendre polynomial of degree k. The splitting 
coefficients with odd k arise from solar internal rotation, while the coefficients with 
even k are caused by departures from spherical symmetry in solar structure, or from effects 
of magnetic helds. Most of the m-averaged power spectra that we £t for this manuscript 
were generated using only the three lowest, odd-/c frequency-splitting coefficients (i.e., oi, 03 , 
and 05 ). In the second step of the process in which we computed the m-averaged spectra, 
we shifted each tesseral and sectoral spectrum by the calculated frequency shift for that m- 
value. In the third step in this procedure, we averaged all of these shifted spectra together 
with the un-shifted zonal spectrum, o(i^), to create the m-averaged spectrum <h™®'''(z/) for 
each degree /. If this m-averaging is carried out in an unweighted manner, we will refer to 
the resulting set of m-aver aged spectra a s “un weighted, m-averaged spectra”. However, as 
we previously described in [Rhodes et al.l (120011) . it is also worth considering average spectra 
computed by combining the 2 / -|- 1 spectra in a weighted manner, using as weights 
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the inverse mean power over the 1500 to 4500 /rHz freqnency range. We will refer to snch a 
set of spectra as “weighted, m-averaged spectra”. 

We have developed two different versions of the cross-correlation method to generate 
rotational splitting coefficients which then are used in the computation of the m-averaged 
spectra. In the hrst of these versions, we cross-correlate the individual spectra over a wide 
range of frequencies such that most or all of the ridges at a given degree are included, while 
in the second version we carry out the cross-correlation over a narrow range of frequencies 
centered about a single ridge at each degree. In the second version, we then repeat these 
narrow-band cross-correlations for all of the successive ridges at a given degree in order to 
build up a set of narrow-band splitting coefficients for that degree. In the wide-band version 
of the cross-correlation code we effectively are computing the averages of the frequency 
splittings over all of the adjacent ridges which are located within the frequency limits of 
the cross correlation (typically from 1800 to 4800 /iHz). The splitting coefficients from the 
wide-band procedure are called the n-averaged splitting coefficients, while those from the 
narrow-band version of our code are called non-n-averaged splitting coefficients. 


2.3. Correction for distortions introduced by latitudinal differential rotation 


For the results that we will be presenting later, we generated a set of n-averaged 
frequency-splitting coefficients by cross-correlating the un-averaged power spectra obtained 
from the 7^1996_61 observing run (cf. Tabled]). The odd-order splitting coefficients (i.e., Oi, 
03 , and Os) that we obtained from this procedure are shown here in the left three panels of 
Figure O These raw frequency-splitting coefficients show large discontin uities in the degre e 
range of 200 < I < 240. Similar discontinuities were hrst noticed by iKorzennild (119901) . 
Subsequently, iRhodes et ahl fjl998al) conhrmed the presence of these jumps in MDI observa¬ 
tions. The exact location of the range of I values where these jumps occur depends primarily 
upon the duration of the observing run from which the power spectra were generated, with 
shorter-duration observing runs showing the jumps at lower degrees. 


WoodardI (119891 ) pointed out that the distortion of high-degree /- and p-mode eigenfunc¬ 


tions caused by a slow, antisymmetric differential rotation can be expressed as a superposi¬ 
tion of the unperturbed eigenfunctions of the same radial order n if the C oriolis forces are 


neglec t ed. Following a dis cussion of Woodard’s suggestion in a preprint of lKorzennik et ah 


(120041) . iReiter et al.l (120031) found that the inclusion of this effect in the calculation of the leak¬ 
age matrices had a very dramatic impact upon the resulting frequency-splitting coefficients. 
Examples of the changes introduced into t he odd splitt i ng co efficients when corrections are 
made for the distortion were presented bv iReiter et al.l (120031) . who showed that at the de- 




















10 


grees below I ~ 200 the splitting coefficients remained almost unchanged, while at the higher 
degrees the jumps were seen to disappear when the mode coupling due to the differential 
rotation of the Sun was taken into account. 


The results that iReiter et ahl (120031) presented were generated using a preliminary 
version of our MPTS method that we have been developing in parallel to the WMLTP 
method that we are presenting in this paper. Because: 1) the MPTS method is extremely 
computationally-intensive; 2 ) it cannot be employed upon power spectra that come from 
observing runs that are as short as only three days in duration due to the low signal-to-noise 
ratios that are inherent in such low-resolution spectra; and 3) we have not yet had the op¬ 
portunity of implementing the changes that we have recently made in our WMLTP method 
into the MPTS code, we have not yet employed the MPTS method to compute entire sets of 
frequency-splitting coefficients. Instead, we corrected the raw splitting coefficients that are 
shown in the left three panels of Figure [2] with a two-step adjustment procedure. First, for 
each of the hve splitting coefficients we £t a least-squares straight line to the degree range 
of 90 to 190 and we also fit a second least-squares straight line to the degree range from 
230 to 400. For the degrees ranging from 200 to 230 the two linear hts for each splitting 
coefficient were simply connected with a third straight line. The difference between that line 
and the extrapolation of the left-hand line was subtracted from the raw coefficient values 
for all of the degrees between 200 and 230. For all degrees above I = 230 the offset em¬ 
ployed for I = 230 was subtracted from each coefficient. Because these initially-corrected 
splitting coefficients showed evidence of systematic variations with increasing degree in the 
three odd-order coefficients (i.e., Oi, 03 , and 05 ), we computed, in the second step of our 
correction procedure, a low-order polynomial ht to each of the three odd-order coefficients 
over the degree range of 488 to 1000. We then subtracted these polynomial fits from the 
partially-corrected, odd-order splitting coefficients and we stopped the correction process 
at this point. This procedure generated the set of corrected odd-order splitting coefficients 
that are shown in the right three panels of Figure IH where it is clear that the jumps have 
been removed. We refer to this set of adjusted frequency-splitting coefficients as our set of 
“corrected, n-averaged” coefficients. 


Using a similar adjustment procedure we also have corrected the set of raw non-n- 
averaged frequency-splitting coefficients that we previously computed using the narrow-band 
version of our cross-correlation method on the un-averaged power spectra obtained from the 
7?.1996_61 observing run. We refer to this set of adjusted frequency-splitting coefficients as 
our set of “corrected, non-n-averaged” coefficients. 
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2.4. Sets of m-averaged power spectra generated from the observing runs 

used in this work 

From the un-averaged power spectra obtained from the observing runs specihed in Ta¬ 
bled] we have generated a total of seven different sets of m-averaged power spectra, which are 
listed in the second column of Table |2] using a naming convention quite similar to that intro¬ 
duced in Table [T] to refer to the individual observing runs. These seven sets of m-averaged 
spectra differed in four key ways: 1) origin of the un-averaged power spectra, 2) whether 
or not the raw, un-averaged spectra were weighted prior to being averaged, 3) whether the 
set of frequency-splitting coefficients was used as computed or as corrected for the effects of 
latitudinal differential rotation and other features that appeared not to be solar in origin, 
and 4) whether or not the set of frequency-splitting coefficients was computed using a narrow 
or a wide frequency range at each degree (i.e., whether those coefficients were computed for 
the individual ridges or were computed in an n-averaged manner). We note that we will not 
present any fits to the m-averaged spectral set iS1996_61 in this work. Rather, this set of 
m-averaged spectra is listed in Table [2] only for the sake of completeness because it was a 
by-product of the cross-correlation process that generated the raw, uncorrected, n-averaged 
splitting coefficients from observing run 7^1996_61. 


3. Problems requiring the use of multiple peaks in the fitting profile 

3.1. Basic considerations 

For the following reasons, high-degree modes cannot be observed as isolated, sharp peaks 
but only as ridges of power. First, the power spectrum computed for a specific target mode 
with degree I and azimuthal order m contains contributions of power from modes with neigh¬ 
boring I and m because the spherical harmonic functions used in the spatial decomposition 
of the observed Dopplergrams are not orthogonal on that part of the Sun we observe (see 
equation ([T])). These unwanted contributions, or spatial leaks, are quantified by the so-called 
leakage matrix (see Sect. 13.21) . Second, with increasing degree the frequency separation of the 
spatial leaks decreases, while the mode linewidth increases with both frequency and degree. 
As a consequence, individual modal peaks blend together to form ridges of power. Typically, 
modes begin to blend into ridges for degrees ranging anywhere from I ^ 20 {n = 29) to 
I ~ 300 (n = 0) depending on the radial order n. Since the amplitudes of the spatial leaks 
are asymmetric with regard to the target mode the central frequency of a ridge is signifi¬ 
cantly offset from the target mode frequency. Therefore, the distribution of power in a ridge 
cannot be simply represented by using just a single symmetrical or asymmetrical function 
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of frequency. Rather, a sum of individual overlapping profiles must be employed the relative 
amplitudes of which are governed by the leakage matrix appropriate to the targeted mode. 
Thus, the correct estimation of the leakage matrix is crucial in the accurate measurement of 
high-degree mode parameters. Moreover, the use of a model prohle consisting of the sum of 
individual prohles allows the htting of low-, medium-, and high-degree modes in like manner. 
In this way systematic errors are avoided which otherwise would be inevitably introduced if 
subsets of mode parameters are to be combined each of which has been generated by using 
a different htting methodology. 


3.2. Leakage matrix 

While the determination of the leakage matrix is straightforward for low- and medium- 
degrees, at high degrees the leakage matrix calculations are greatly complicated by the 
necessity to take into account (1) the horizontal component velocity, (2) the distortion of 
the eigenfunctions by the solar differential rotation, and (3) instrumental effects that cause 
image distortion and smearing. We will address these issues in the following paragraphs. 


3.2.1. Radial and horizontal component 

The Fourier transform Oi^m of the time series of spherical harmonic amplitudes of a 
target mode with degree I and azimuthal order m can be written as a sum over the Fourier 
transform dn>,i>,m> of the time series of the solar oscillation modes given by viz. 


Ol,m ^ ^ ,m'On>p^rri '1 

n' ,1' ,m' 


(5) 


where is the leakage matrix. As iKorzennik et ahl fl2004j) have shown, the leakage 

matrix can be written as 


^ m I,m;l',m'^ 

where is the part coming from the radial displacement, is the part coming 

from the horizontal displacement, and ^ is the ratio of the horizontal displacement to the 
radial displacement. It should be noted that both and independent of 

the radial order n' of the mode. Using the normalization of the displacement eigenfunction 





~ 13 - 


components given bv iGoughI fjl993l) . the displacement component ratio is given by 


{n,l) ^ 

“id?wi 


(7) 


Here, r* is the radial location of observation, and and are, respectively, the horizontal 
and radial components of the displacement eigenfunction of the mode (n, /), given by 


^{r,e,(p) = 




(h). 


+ ^ L,l 


d 


d 


'(p 


sin 9 d(j) 




( 8 ) 


where (r, 0, 0) are spherical polar coordinates with r being the distance to the center, 9 being 
the co-latitude, 0 being the longitude, e^, e^, e<^ are, respectively, the unit vectors in the 
r, 9, (j) directions, Yim{9,(i)) is a spherical harmon ic of degree I and azimuthal order m, and 
= l{l + 1). As IChristensen-DalsgaardI (120031) has shown, the displacement component 
ratio can be written as 


{n,l) _ GMqL _ 


Cl •' = 




where 


_ [GMqL 


1 ^ 0,1 

l^n,l 

1/2 


(9) 


( 10 ) 


is the frequency of the /-mode of degree I in the asymptotic high-degree limit (lGoughlll980l) . 
Un,! is the average frequency for the multiplet (n,/), and Mq and Rq are, respectively, 
the mass and the radius of the Sun. For p-modes equation ([9]) implies that 0 < q 


(n,l) 


t ^ 


< 


1 because, for hxed /, the /-mode frequency is smaller than any p-mode frequency. It 
should be noted that as dehned in equation ((Tj) is not only equal to the ratio of the 
displacement eigenfunction components but is also equal to the ratio of the horizontal and 
vertica l components of the v elocity eigenfunctions for the mode. For low- and medium-degree 
modes Rhodes et ah ( 2001 1 have shown tha t the observed values of closely match the 


theoretical prediction given in equation (|9]). iRabello-Soares et al.l (l200ll ) arrived at a similar 
conclusion. For high-degree modes the agreement between the measured and theoretica l 


horizontal-to-vertical displacement ratio has been demonstrated bv ISchou fc BogartI (119981) . 


When power spectra are to be htted rather than Fourier spectra we have to compute 
the leakage matrix, relevant to power spectra, which is given by 


(ps) _ /^2 


'' ,7n' 


,m' * 


( 11 ) 


Moreover, for the htting of m-averaged power spectra we have to compute the leakage matrix, 
which measures, for a given ridge of radial order n, the contribution of a mode of 
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given V in the power spectrnm calculated for a mode of given 1. To do so, we need to take 
the sum of the squares of all the m-leaks in equation ([6]). Using equation flTTl) we get 


i v 

^(mavg) _ ^(ps) 

/ j / j ‘ 

m=—l m'=—V 


( 12 ) 


In practice the leaks fall off rather rapidly with increasing |m — m'| and increasing 

\l — l'\. Hence, the sums in equation flT^ only have to be evaluated for a limited range of 
both |/ — l'\ and \m — m'\. 


3.2.2. Distortion by the solar differential rotation 


One of the conspicuous effects of solar rotation is the well-known splitting of the oscilla¬ 
tion frequencies, that is the dependency of the oscillation frequencies on the azimuthal order 
m (cf. equation (|1])). Similarly, the modal eigenfunctions of the solar oscillations depart from 
their custom arily assumed sp herical harmonic form (cf. equation (|8])) as a result of solar 
rotation. As IWoodardI (119891) has shown, the distortion of high-degree mode eigenfunctions 
by a slow, axisymmetric differential rotation can be expressed as a superposition of the un¬ 
perturbed eigenfunctions of the same radial order n, if Coriolis forces are neglected. He also 
has shown, that the perturbed leakage matrix can be expanded in terms of the unperturbed 
leakage matrix as 

(13) 


where 


-1? 


7z/,z// = <^ 271 

0 


cos [px + Z{x)] dx for I' — I" even 
for I' - V odd 


(14) 


p={i'-n/2 


Z{x) = 


xf, m 


4 dn/dl\ii 

9 

m 


(i ?2 + sin ^ -sin (2x) 


(15) 

(16) 

(17) 


In equation ffTBD dv/dl\ii denotes the derivative of frequency v with respect to degree I 
evaluated at degree I'. It is assumed that I can be treated as a continuous variable. Actually, 
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dv/dl is calculated by taking the derivative of a smooth function htted to the march of 
V versus I along a ridge of given radial order n (cf. Section 


The ^-coefficients in 

equation fITB]) result from a parametrization of the angula r velocity VL(9) of t h e sur face 
differential rotation as a function of co-latitude 9. Following ISnodgrass fc Ulrich! fjl990f) we 
have 

n{e) = 27r(5o + B 2 cos^ 9+ Bi cos^ 9), (18) 


where 

Bq = 473.0 nHz, B 2 = —77.0 nHz, i ?4 = —57.5 nHz 
for rotation of Doppler features on the solar surface. 


(19) 


The dependence of the perturbed leakage matrix flT^ on the ^-coefficients involved in 
the rotational model ffTSj) implies the following problem. On the one hand, the radial varia¬ 
tion in the latitudinal differential rotation profile and, hence, the surface differential rotation 
can be measured through a rotational inversion of the frequency-splitting coefficients 
as given in equation (jTj). For the measurement of the a perturbed leakage matrix flTSD 
must be specified and input into the peak-bagging code. On the other hand, equation ffT6|) 
demonstrates that the perturbed leakage matrix itself and, hence, also the frequency-splitting 
coefficients, depend upon the il-coefficients which were used to parametrize the solar differ¬ 
ential rotational profile in the surface layers. This mutual dependence can only be resolved 
with some kind of fixed-point iteration. Nevertheless, as we discussed earlier in Section 12.31 
the importance of including these e ffects into the calculation of the leakage matrices was 
demonstrated by Reiter et ah fj2003 h who showed that their inclusion removed the disconti¬ 
nuities that were otherwise present in the non-n-averaged splitting coefficients for the n = 2 
ridge. 


3.2.3. Instrumental effects 


For the determination of high-degree mode parameters not only the leakage matrix must 
be known but also the instrumental characteristi cs must be very well understood and very 
precisely measured fcf. lRabello-Soares et al.ll200lh . The instrumental effects to be considered 
in the analysis include plate scale error, image distortion, width and spatial non-uniformity 
of the instrumental point spread function (PSF), image orientation (P-angle), and the finite 
pixel size of the detector. Similarly, one has to consider errors in the P-angle and Po-angle 
caused by errors in the ass umed orienta ti on of the solar rotation axis. Using data obtained 
with the MDI instrument, iReiter et al.l (120031) have shown that the inclusion of both the 
plate scale error and the image distortion has a rather strong impact upon the measured 
splitting coefficients for degrees I > 200. 
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We have found it convenient to calculate the leakage matrix by constructing simulated 
images corresponding to the line-of-sight contribution of each component of a single spherical 
harmonic mode and then decomposing that image into spherical harmonic coefficients using 
exactly the same numerical decomposition pipeline employed to process the observations. 
This approach has the advantage that some of the above-described instrumental effects can 
easily be included in the leakage matrix calculation. It should be noted, however, that the 
inclusion of instrumental effects in the effective leakage matrix is, in general, not equivalent 
to taking into account those effects in the pipeline used to process the observations. 


For a th orough discussion of the inst rumental effects of the MDI instrument we refer 
the reader to Korzennik et ah ( 2004 . 2008 ). 


4. The windowed, multiple-peak, averaged-spectrum method 


The Windowed, MuLTiple-Peak, averaged-spectrum method, which we will refer to in 
the following as either Method 2 or the WMLTP method is an advancement of our Method 1 
that is briefly described in Appendix Several steps were involved in coming to the design 
of the new method. As compared to Method 1 the improvements incorporated in Method 2 
include (1) the replacement of a single symmetric prohle with a sum of asymmetric prohles 
representing the target peak as well as the peaks of the neighboring /-leaks, (2) a sum of 
asymmetric prohles representing the peaks of the n-leaks neighboring the target peak, (3) 
the temporal side-lobe peaks, and (4) an approximation to the m-averaged leakage matrix. 
Generally speaking, an /-leak is a spectral peak of the same radial order as that of the target 
mode but whose degree is different from that of the target mode, while a n-leak is a mode of 
arbitrary degree whose radial order is different from that of the target mode. We note that 
all of the results that we will present in this paper were obtained with the latest version of 
the Method 2 code, which we refer to as the rev6 version of this code. 


By the time that we began to develop our WMLTP htting method, the observations 
of iDuvall et ahl fll993l) had indicated that the peaks in the solar oscillation power spectra 
showed varying amounts of asymmetry. In particular, Duvall and his collaborators noted 
that the velocity and intensity power spectra revealed an opposite sense of asymmetry. 
After these observations appeared, we (IRhodes et al.lll9971) modihed our Method 1 htting 
method in an ad hoc manner by employing an as ymmetric pro h k tha t consisted of two 


Lorentizan half-prohles that had differing widths. iRhodes et ahl fjl997l) used that ad hoc 


prohle to demonstrate that the observed asymmetries in the peaks in MDI velocity power 
spectra shifted the htted frequencies by a substantial amount in the fre quency range where 
structural inversions are most sensitive to the observed frequencies. The iDuvall et ahl fll993l) 
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observ ations were later confirmed by iNigam et al.l fjl998l) . and then iNigam fc Kosovichev 
fll998l) introduced a theoretical prohle that accounted for the observed asymmetry in a self- 
consistent manner. While, in principle, we could have replaced the split-Lorentzian prohle in 
Method 1 with the INigam fc Kosovichevi fjl998l) prohle and continued to use that modihed 
method, the problems that we described above in Section [3] caused us to abandon that 
method as well as our use of the ad hoc split Lorentzian profi l e. In stead, we continued 
to develop our WMLTP method using the INigam fc Kosovichevi (119981) prohle. As we will 
note below, the INigam fc Kosovichevi (119981) prohle turns into a symmetric Lorentzian prohle 
when its asymmetry parameter is set equal to zero; hence, its adoption allows us to ht power 
spectra using both asymmetric and symmetric prohles. 


4.1. Implementation of m-averaged leakage matrices 


As we have already mentioned in Section 13.2.31 we are calculating the leakage matrix 
(cf. Section 13.21) by using the very same numerical decomposition pipeline that is employed 
to process the observations. This approach provides initially the leakage matrix elements 
froui which we then compute the m-averaged leakage matrix by means of 

equations (ITT]) and (IT^ . 

As is demonstrated here in Figure [3] the m-averaged leakage matrix, can be 

approximated, for given radial order n and degree I, by a Gaussian prohle, viz. 


a 


(mavg) 


n\l^V 


= exp (A/ - > 0, (20) 


where is a parameter related to the width of the leakage matrix, Xm’^'^ is a parameter 
that accounts for the ohset of the leakage matrix due to the horizontal component of the 
modal velocity eigenfunction of the Sun (cf. Section [3. 2. ip . and Al = I' — I is the distance of 
the spatial leak located at degree I' from the target mode. In both panels of Figure [3] the m- 
averaged leakage matrix is marked by diamonds, while the htted Gaussian prohle, as given in 
equation (l20p . is represented by the full line. Overall, in both panels the march of the leakage 
matrix is well described by the Gaussian prohle. The largest deviations are just in the ten 
percent range, and moreover do occur at locations at which the amplitude of the leakage 
matrix is small. Hence, we believe that those deviations can safely be neglected. In Figure [3] 
we also show in both panels the leakage matrix that results if the distortion introduced 
by the solar latitudinal diherential rotation is taken into account (cf. Section [3.2.21) . The 
corrected leakage matrix is marked in both panels by the triangles, and the corresponding 
htted Gaussian prohle, as given by equation (l20l) . is represented by the dashed line. We note 
that for the (n, 1) = (2,1000) mode the width of the corrected leakage matrix greatly exceeds 
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the width of the uncorrected leakage matrix, while for the (0,400) mode the corrected leakage 
matrix is only slightly wider than the uncorrected leakage matrix. While the Gaussian prohle 
represents the march of the corrected leakage matrix quite good for the (0,400) mode, in the 
corrected leakage matrix for the (2,1000) mode strange asymmetric “shoulders” do appear 
at about Al = —5 and A/ = +4, respectively, that are not well represented by a simple 
Gaussian prohle. We hope to re-visit the issue of the suitability of the Gaussian prohle in a 
later version of the WMLTP method. 


The Gaussian approximation, as dehned in equation fl2U]) . to the m-averaged leakage ma¬ 
trix is implemented into the WMLTP method by means of the following approach. For given 
degree I and radial order n, we hrst compute for a set of discrete values of the displacement 
component ratio (cf. Section [3.2.ip . 


= j/10, j = 0,..., 10, 


( 21 ) 


the m-averaged leakage matrix that has been corrected for the latitudinal diherential 

rotation as described in Section 13.2.21 Next, we ht the Gaussian prohle, as given in equa¬ 
tion (120]) . to the resulting leakage matrices to get, for each of the values j = 0,..., 10, 
the values and of the ht parameters and Xm'^\ respectively, that determine 
the Gaussian prohle, as dehned in equation fl20l) . Because the three parameter s Xm'^\ 

and are not independent from one another, we follow Rhodes et al. (2001) and expand 
both and 4"’^^ in terms of Xm'^\ viz. 


(n,Z) ^ ^(n,Z) ^(n,Z) ( 



(22) 


(23) 


where and denote the respective expansion coefhcients. 

We note that it would be mathematically equivalent to expand both aGl) Xm’^'^ in terms 
of In the hnal step of our approach we ht the expansion (|22|) to the knots > 

j = 0,..., 10, and the expansion (l23l) to the knots j, j = 0,..., 10, to get the set 

of expansion coefhcients and Typical hts of both 

and versus Xm'^^ are shown here in Figure IH for the modes (n, /) = (11,10), (0,400), and 
(2,1000), respectively. The hts clearly demonstrate that the expansions given in equations 
(]22|) and (l23l) . respectively, are reasonable approximations to the variation of both and 
with respect to Xm’^\ Moreover, from the left panel in Figure [Hit becomes evident that 
q,()t.,z) decreases with increasing degree /. Hence, according to equation (120|) the width of the 
m-averaged leakage matrix increases with increasing degree /. 
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Once the set of expansion coefficients and has been 

computed for the mode {n, 1), the corresponding m-averaged leakage matrix Anj, as given in 
equation (120|) . can easily be determined by performing the following steps. First, we calculate 
the value of from equation ([9]). Second, using this value of we solve equation fIMl) 
for Finally, we insert this value of into equation fl22ll to get the value of 


4.2. Theoretical fitting profile 


In the WMLTP method the following fitting profile is used to represent an oscillation 
peak of given degree I and radial order n: 




(24) 


Here, u is frequency, W(i^) is the Fourier spectrum of the temporal window function of the 
observational time series, “(g)” denotes the convolution operator, and AAn,i{^,p) is defined 
by 


M 


N 


MnA^,p) = ^(A^,i,B^,i,XAi)e(Al)A^A^l,A^’^\xi:’‘^) + '^^(A,B„x,), 


Al=-M 


i=l 


where 


targeted peak and /-leaks 

+ a + bi' + ci'‘^, 

V* 

background 

1 4 - 


n-leaks 


^{A,B,x) = A 


/^ma 


At=l 


- "nJ + 2j (^0 


(25) 

(26) 
(27) 
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“ dl^ 
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n,l I 1 + ^ 
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(28) 


(29) 
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Here, ^ is the asymmetric profile of iNigam fc Kosovichevi fjl998l) . 0 is an empirical adjust¬ 
ment to the amplitudes of the /-leaks, ^ is an approximation to the m-averaged leakage 
matrix that is described in Section H7T1 2M is the number of the /-leaks, N is the number of 
the n-leaks, Hj, Bi, z/j, and Wt are, respectively, the amplitude, the line asymmetry param¬ 
eter, the frequency, and width of the i-th n-leak, i = 1,... ,N, and a, b, c are parameters 
describing the background noise. We note that in equation (|2^ it has been presumed that 
each of the 2M +1 modal peaks can be characterized by the same line asymmetry parameter 

Bn,l- 


The total number of /-leaks, 2M, included in the model prohle, as given in equation fl25|) . 
depends on the width of the leakage matrix, An^i, and is typically in the range from 10 to 
60. As we have noted in Section ITTl particularly for higher degrees, /, the leakage matrices 
are rather wide which is mainly due to the distortion of the oscillation eigenfunctions caused 
by the latitudinal differential rotation. Therefore, in equations (127|) and (128|) Taylor series 
expansions of fairly high order have to be employed to adequately describe the variation of 
amplitude, frequency, and linewidth with degree /. In the current version of the WMLTP 
code /imax = 8 is used. For the determination of the total number of n-leaks, N, to be 
included in the model prohle given by equation (I25p . we refer the reader to Section 031 In 
our applications we found that N can vary from 0 up to about 23, depending on both degree 
/ and radial order n. 


In the htting prohle dehned in equations fl2T|) through fl29]) a total of 4A^ -|- 8 htting 
parameters are involved, viz. the mode amplitude A„the mode frequency the mode 
linewidth Wn,i, the background noise parameters a, b, c, the line asymmetry parameter Bn^i, 
the ohset, Xm'^\ of the m-averaged leakage matrix due to the horizontal component of the 
modal velocity eigenfunction of the Sun (cf. Sect. 14.ip . and the parameters Aj, i/j, tCj, Hj, 
i = 1,..., A^, representing, respectively, amplitude, frequency, linewidth, and line asymmetry 
of the N n-leaks. These htting parameters can be lumped together in the htting vector 


T 

p = Wn,h Bn,h ^ 0 Ai, Ui, Wi, Hi, ... , A^r, Wn, Bn) . (30) 

We did not include in this htting vector p the parameter because it depends on 
by virtue of equation ([22]). We also did not include the Taylor series expansion coefficients 


1 d^A 1 d^w ^ 

AIF' W w1»’ = 


(31) 


for the amplitude, frequency, and linewidth, respectively. Rather, these coefficients are taken, 
for given degree / and radial order n, from tables containing initial estimates or so-called 
seeds that are derived from previously computed modal parameters, a n d are improved upon 


in a hxed-point iteration similar to that described in [Rhodes et ahl fl200ll) . As a result. 
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the compilation of such seed tables is an important pre-processing step for the use of the 
WMLTP method. 

In order to investigate the performance of our implemented hxed-point iteration we have 
done the following analysis in which we have restricted ourselves to the investigation of the 
frequency, u, and the derivative thereof with respect to degree, Au/Al. For convenience, 
we did not consider either the linewidth or the amplitude and the derivatives thereof with 
respect to degree. In the hrst step of this analysis we used our seed table for the epoch of 2010 
which provided both the seed frequencies, and the derivative thereof with respect 

to degree, to £t the m-averaged spectral set iS2001_90. From this htting 

run we obtained the frequencies, ■ We chose the m-averaged spectral set 52001 _90 for 
this test because the 7?.2001_90 observing run corresponded to the maximum phase of Solar 
Cycle 23. Therefore, due to the well-known shifts in the /- and p-mode frequencies that 
accompany changes in solar activity, the htted frequencies, and the frequency derivatives, 
that would result from the htting of these power spectra, would be expected to exhibit the 
largest possible differences from the frequencies and derivatives in our 2010 seed table since 
the latter were generated from the 7?.2010_66 observing run when the mean level of solar 
activity was much lower. 

Next, we fitted to the frequencies, on a ridge-by-ridge basis, a smooth function 

of degree, to get new seeds and respectively, by setting 

(32) 

(33) 

i 

where the superscript “(n)” denotes the ridge of radial order n. For the details of the 
construction of the function, we refer the reader to the end of Section [6.21 In the 

second step we used this new set of seeds to £t the m-averaged spectral set 52001 _90 once 
again to get the frequencies , from which we generated yet another set of seeds 
and {Au/ respectively, in the same manner as we have done in the previous step. 
These new seeds were then used to £t, in the third and hnal step, the m-averaged spectral 
set 52001_90 a third time to get the htted frequencies, ■ For further explanation we 
have illustrated the sequence of individual steps performed in our analysis in Figured 

The results of our analysis are summarized in Table [3l As can be seen from columns 4 
and 5, the two-fold update of the seed table resulted in a signihcant reduction in both the 
average and the standard deviation of the diherences of the seed values of Au/Al as well as 
of the scaled diherences of the htted frequencies. As expected, the rather large changes in 


^^seed,2) ^ 

seed,2) ^ Qp{n)(^l-^ 


dl 
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the scaled differences of the htted frequencies that accompanied the hrst update of our 2010 
seed table were reflective of the very large difference in the mean levels of activity between 
2001 and 2010. The fact that we could obtain a self-consistent seed table for 2001 with only 
two iterations in spite of such a large difference in the level of activity in the two epochs 
justihes our approach of replacing the Taylor series expansions, as given in equation ([3T]), 
with such self-consistent seed values. 


4.3. Selection of fitting box widths 


The width of the selected htting box is crucial for the successful determination of the 
htting vector p dehned in equation (l30D . We have found it useful to construct the htting 
box for the mode {n, 1) as follows: 


low _ seed _ {n,l) 

^n,l — ^n,l ^ 


^^n,l 

An 


left 


— ,,seed I (n,l) 
^n,l — ^n,l + ^ 


^^n,l 

An 


right 


(34) 


Here, is the seed frequency of the targeted mode (n,/), and are the lower and 
upper boundary, respectively, of the htting box, An„ ;/An|ieft and Az/„;/An|right denote the 
variation of the mode frequency with respect to the radial order n at the left (i.e., lower 
frequency) and right (i.e., higher frequency) side, respectively, and is given by 
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(35) 


where and 4"^^ are predetermined parameters depending on the radial order n. 

The quantities Anrr^;/An|ieft and Az/„^)/An|right in equation fl34)) are approximated as follows: 
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for n > 0, 
for n = 0, 


(36) 


where n®®®*^, and n®®^*^; are taken from a seed table. As to the choice of the parameters 

4”^ in equation fl3^ it must be kept in mind that any selected htting box 

must fulhll at least two requirements. First, the htting box must be sufficiently wide so that 
both the mode prohle and the background power are well sampled. This requirement becomes 
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an issue if spectra are to be fitted that are derived from an observing run the duration of 
which is only a few days. In this case it can happen that the number of £t parameters 
included in the htting vector p exceeds the number of frequency bins constituting the htting 
box. Second, the fitting box must not be unduly wide in order to save computing time 
which non-linearly increases with the number of frequency bins comprising the htting box. 
In practice, we determine the parameters and for a given ridge of radial 

order, n, by a trial-and-error method and select those values of them which give rise to the 
least scatter of frequency along the ridge. 

As can be seen from equations flM|) and fl35|) . aside from variations of both AAn|left 
and An„^i/An|right with degree /, the width of the fitting box is constant for I < and 
I > and varies linearly with degree I for < I < We note here that in general 
An„i/An|ieft 7 ^ Az/„ ^/Anlright- Therefore, according to equation flM|) the fitting box is 
generally not symmetric with respect to the seed frequency 

For a selected set of radial orders, n, we show in Figure [6] the half of the width, Wb, 
of the respective htting boxes measured in terms of the average, (Au/An), of Az/„ i/An|ieft 
and An^^i/Anlrighf. For most of the cases shown the width of the htting box increases with 
increasing degree, while for n = 6 the width is practically constant, and for n = 8 and n = 29 
the width decreases with degree. We also note that for the higher-order ridges the n-leaks 
(n ± 1, /) are included in the htting box. For the n = 20 ridge the htting box is getting so 
wide that even the n-leaks {n ± 2, 1) are encompassed. The choice of such wide htting boxes 
is expressive of the fact that, in general, the scatter of the fitted mode parameters along 
a given ridge decreases with increasing width of the htting boxes. Therefore, it would be 
desirable to ht all n-values simultaneously that are present in an m-averaged spectrum of 
given degree. However, such approach is not feasible on practical grounds. 


4.4. Determination of the number of n-leaks to be included in the fitting 

model profile 

The number of n-leaks, N, to be included in the htting model prohle dehned by equa¬ 
tions (|2T|) through ([29]), cannot be determined from hrst principles. Therefore, a range 
1^ — ^o| A Grange, ” ''^ol A ''Grange, of I and u values is choseu around the targeted mode 
(no,/o)) cind every mode (n,/) within this range that (1) has a frequency within the fre¬ 
quency range of the htting box (cf. Section IT73|) . (2) dihers from (no,/o )5 and (3) is not 
identical to any of the /-leaks, is included as an n-leak in equation fl25|) . In the current 
version of the WMLTP code we are using /range = 4, rirange = 3 for isolated modal peaks, 
(uq, /q), and /range = 4, Urange = 2 Otherwise. Usually, many of the n-leaks determined in this 
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manner turn out to be statistically insignificant for getting a reasonably good fit. Because 
generally the n-leaks included in the fitting profile strongly affect the fitted parameters of 
the targeted mode {rioJo), it is crucial to discard the statistically insignificant n-leaks while 
keeping the statistically significant ones. Therefore, the statistical significance of an n-leak 
is tested by virtue of the so-called R-test flFriederu 119831) . This test is skipped, however, for 
“true” n-leaks (no ± 1, 1) and (no ±2,1), respectively, with (no, /o) being the targeted mode. 
Because the R-test does not always work satisfactorily, we have additionally implemented 
into the WMLTP code heuristic criteria for discarding an n-leak. For example, an n-leak is 
discarded if its amplitude and/or width is outside a given range or if it overlaps too much 
with another n-leak or with the targeted peak itself. Overall, the determination of the num¬ 
ber of n-leaks, N, is a rather time consuming undertaking which significantly increases the 
compute time of the WMLTP code. 


4.5. Implementation of numerical scaling, adjustments, and options 

The vector p defined in equation fl30l) is determined by fitting, in the least-squares 
sense, the model profile given in equations fl2T|) through fl29l) to the m-averaged spectrum 
of given degree I in a fitting box (cf. Section 14.31) centered about the target peak of radial 
order n. The confidence interval on each fit parameter pi E p is calculated as described 
in Appendix In order to make the solution of this least-squares problem as robust as 
possible diverse provisions are made. First, all fit parameters are scaled such as to make 
their values on the order of unity. This scaling greatly improves the numerical stability 
of the solution of the least-squares problem. Second, the fit parameters representing the 
target peak (viz. h'n,^ Wn,i, B^^i) as well as the fit parameters representing the n-leaks 
(viz. Ai, Pi, Wi, Bi, i = 1,...,A^) are subject to bounds, i.e., they are constrained to he 
within prescribed intervals. This is important in order to avoid unphysical values of any of 
these fit parameters, for example, negative amplitudes, frequencies too much off from their 
seed values, linewidths much smaller than the spectral resolution, or linewidths on the order 
of the width of the fitting box. Also, the magnitude of the line asymmetry parameter of 
both the target peak and the n-leak peaks is constrained to values less than or equal to 0.3. 
Third, the fitting vector p is determined in several passes. In Table H] it is shown which 
parameters are fitted in the individual passes. Any parameter not listed in a given pass as 
a fitted parameter in Table IU but which is included in the fitting vector p, is kept fixed 
either to its seed value or else to its value obtained in a previous pass. We note that the 
permutations given in Table 0] were chosen on purely heuristic grounds and are, therefore, 
somewhat arbitrary. Fourth, a flag can be set in the WMLTP code which enforces the use of 
the symmetric Lorentzian profile. This is accomplished by setting R = 0 in equation ([26]). 
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Otherwise, the full asymmetric prohle is used, in which case the asymmetry parameter, B, 
in equation fl26|) is employed as a htting parameter in the least-squares problem. 


For the solution of the nonlinear constrained least-squares problem which results for 
the determination of the htting vector p we use the FORTRAN code NLPQL which is a 
special implementation of a so-called sequential quadratic programming metho d, and which 
is gen erally designed for solving constrained non-linear optimization problems fjSchittkowski 
19861 ). Sequential quadratic programming is one of the most robust algorithms for the 


solution of non-linear continuous optimization problems. The method is based on solving a 
series of sub-problems designed to minimize a quadratic approximation to the Lagrangian 
function subject to a linearization of the constraints. 


In the case of either high frequencies and/or high degrees when the individual modal 
peaks can no longer be resolved but rather blend together to form ridges of power in the 
m-averaged spectra, the determination of the htting vector p given in equation fl30|) becomes 
an ill-dehned least-squares problem. This is because in this case any change of the parameter 
does not result in a signihcant change of the htting prohle dehned by equations 
through (12^ because which is related to the width of the leakage matrix, turns out to 

be practically independent of as is demonstrated here in the left panel of Figure HI As 
a result, Xm’^'^ can no longer be used as a ht parameter but rather must be kept hxed to its 
seed value. 


To make the WMLTP method more hexible in practical applications we have imple¬ 
mented the following features into the code. First, if frequency splitting coefhcients are 
available for a mode {n,l), the m-averaged spectrum can be calculated from the respective 
zonal, sectoral, and tesseral spectra in a weighted or unweighted fashion within the htting 
box of that mode making superhuous the calculation of the m-averaged spectrum in the 
pre-processing step as described in Section [221 In other words, this option makes it possible 
to either use splitting coefhcients or else an m-averaged spectrum as input to the WMLTP 
code. This feature is particularly useful if non-n-averaged splitting coefhcients are available. 
Second, for a mode (n, 1) the expansion coefhcients which are required to calculate the corre¬ 
sponding m-averaged leakage matrix parameters and via equations fl22]) and fl22]) . 
can either be read from a hie generated in a previous calculation or else can be calculated as 
part of the htting procedure itself. This way a considerable amount of compute time can be 
saved. Third, if the spectral peak of a mode (n, /) is well separated from the I- and n-leak 
peaks, the ohset parameter of the leakage matrix, Xm'^\ can be invoked as a ht parameter. 
By then conv e rting the htted values of Xm’’’'^ into values of by means of equation 


Rhodes et ahl fj200l[) were able to demonstrate, using a previous version of the WMLTP code. 


that the such measured values of match theoretical predictions. Broadly speaking, such 
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approach is possible for modes with Wn,i ^ 0.8 {Au/ Al)n,i, where Wn,i is the mode linewidth 
and {Au/Al)nj is an approximation to the derivative of the mode frequency with respect to 
degree. 


4.6. Enforcement of symmetrical line profile for the fitting of psendo-modes 


We initially employed the asymmetrical prohle for all of the cases along each ridge. 
However, for frequencies u > 5000 /rHz we found unexpected excursions in the line asymme¬ 
try parameter, B, an example of which is shown here for the n = 6 ridge in the left panel 
of Figure [71 Because such behavior of the line asymmetry parameter, B, is hard to explain 
on physical grounds, we rather presume that the asymmetrical prohle becomes invalid for 
modes with frequencies close to or greater than the acoustic cut-off frequency. This pre- 
sumption is substantiated b y the fact that for those frequencies the assumptions made by 
Nigam fc KosovichevI fll998h in the derivation of their asymmetrical prohle are not fulhlled. 
This is because the high-frequency peaks are not normal modes but rather so-called pseudo¬ 
modes caused by interference between the waves coming directly from the excitation source 
and waves refracted in the i nterior. Hence, the pseudo-modes are not osci llations that are 
trapped in the solar interior flNigam et al.lll998l:lNigam 8z Kosovichevlll998h . Moreover, the 
pseudo-mode spectral peaks are essentially symmetric. Asymptotically, at high frequencies, 
their prohle is described by a sin^ function. Unfortunately, there is no simple htting for¬ 
mula describing both normal modes and pseudo-modes. Because we need to leave the rather 
complicated issue of the proper prohle to be employed as a subject of further investigation, 
we rather decided to resort to the workaround to simply force the htting code to switch 
from the asymmetrical to the symmetrical prohle at a prescribed frequency, z/^, along a given 
ridge. This switchover frequency is chosen to coincide with the hrst zero-crossing of the line 
asymmetry parameter, B, that occurs for frequencies u > 4600 /iHz along a given ridge. For 
example, for the n = 6 ridge the switchover frequency is = 4815 /iHz, as is shown here 
in the right panel of Figure [71 Expectedly, a sharp kink is introduced in the run of B at 
the switchover frequency. Of course, a more physically motivated htting prohle would give 
rise to a gradual transition between the asymmetrical and the symmetrical prohle. We point 
out that the switchover frequency Us depends not only on the radial order n but also on the 
epoch of the observation, because the mode parameters are ahected by the mean level of 
solar activity during each observing run. 


At a glance, the tiny error bars as shown in the left panel of Figure [3 seem to indicate that 
the excursions of the line asymmetry parameter, B, to both positive and negative values for 
frequencies z/ > 5000 /rHz are a statistically signihcant ehect. This is not the case, however. 
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Rather this “signihca nt” asymmetry in the f itted line prohles is a further argument that the 
asymmetric prohle of iNigam fc Kosovichevi fjl998l) is invalid for the htting of pseudo-modes. 
Namely, if it were valid, the resulting line asymmetry should be statistically compatible with 
R = 0 for frequencies v > 5000 pHz because it is known that the pseudo-mode peaks are 
essentially symmetric. 


The use of the symmetrical prohle for frequencies above the switchover frequency, 
along a given ridge not only resolves the issue of the unexpected excursions in the line 
asymmetry parameter, R, but also results in a substantial reduction in the frequency scatter 
observed in the high-frequency portion of each ridge. For the study of the variation in the 
frequency scatter along a given ridge we found it useful to evaluate the point-to-point scatter, 
27, as dehned by equation flBip . not for the frequency v itself but rather for the numerical 
derivative of the frequency with respect to degree, Av/Al. An example of the reduction 
in the scatter of Au/Al for frequencies above the switchover frequency is shown here for 
the n = 4 ridge in Figure [H For this ridge the switchover frequency is Ug = 4700 yuHz. 
In the left panel of Figure |8] the asymmetrical prohle has been used for all of the cases 
along this ridge, while in the right panel the asymmetrical prohle has been replaced with 
the symmetrical prohle for frequencies u > Ug, i.e., for degrees I > 562. In order to evaluate 
the obvious reduction in the frequency scatter quantitatively, we computed 27 for the high- 
frequency portion of the n = 4 ridge by setting /min = 562, /max = 1000, and v = Av/Al in 
equation flBip . When we did so, we found after eliminating some obvious gross outliers that 
27 = 0.639 for the asymmetrical prohle, but only 0.214 for the symmetrical prohle. Hence, 
by simply changing over from the asymmetrical to the symmetrical prohle at = z/^, we were 
able to reduce the scatter in Au/Al by a factor of about 3. We also saw similar reductions in 
the scatter in the high-frequency cases of Av/Al when we employed the symmetrical prohle 
in place of the asymmetrical prohle for the high-frequency portions of the /-mode and the 
other p-mode ridges. 


4.7. Determination of the background portion of the theoretical model profile 

A power spectrum computed from a time series of Dopplergrams of length T has a 
spectral resolution of Av = 1/T. If the Dopplergrams are observed at a cadence of At the 
spectrum covers a frequency range from zero up to the Nyquist frequency v^y = 1/(2At), 
which is 8333 /iHz for the MDI instrument, as we have already pointed out in Section 12.11 
On the other hand, the maximum width, IFmax, of the htting boxes used for htting the modal 
peaks in the set of spectra obtained from an observing run is on the order of a few hundred 
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/iHz at most. Hence, for T being greater than a day or so, we have 


Au < H4iax < i^Ny 


max 


(37) 


The inequality H4iax ^ in the above equation flTTI) implies that we can safely approxi¬ 
mate the background noise present in the measured power spectra as a quadratic function 
of frequency within a selected fitting box. We note, however, that generally a linear back¬ 
ground model, i.e., c = 0 in equation fl2^ . is an adequate choice. We also note that the 
background portion of the theoretical model profile, as given in equations fl24|) through ([29]), 
is not a reliable estimate of the actual frequency variation of the background noise within 
the selected htting box around the targeted peak. This is evident from the fact that the 
£t parameters a, b, and c in equation fl2^ . which describe the background portion of the 
model profile, can change significantly from one degree to the next along a given ridge. This 
scatter in the background portion of the model profile is most likely an artificial effect, and 
unfortunately translates into a similar scatter in other fitting parameters included in the 
model profile, such as frequency and linewidth, for example. To mitigate this scatter, we 
devised the following heuristic approach. First, we fit a straight line to the spectral power 
in the troughs of the spectrum in a frequency range that is twice as large than the fitting 
box centered around the targeted peak to estimate the slope of the background noise. This 
is demonstrated here in Figure [9] where we show, for the modes {n, 1) = (24, 8), and (24, 56), 
respectively, the measured power spectrum in black with the such fitted straight line over¬ 
laid in red. While the troughs in a spectrum are not indicative of the background noise, 
but rather correspond to the intersection of the wings of the peaks, we believe that they 
are indicative of the slope of the background. Second, if be denotes the slope of the such 
determined straight line and Abg the uncertainty thereof, we constrain the slope b of the 
background term in equation fl25]) by 


(38) 


be — 3 Abe A b < be + 3 Abe, 


where the term 3 Abg allows for uncertainties in the value of be- Third, we constrain the 
background term in equation fl25|) by 


a + bu + cu'^>0 for all frequencies u within the htting box. 


(39) 


just to ensure that the estimated background noise is positive throughout the entire fitting 
box. 
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5. Sensitivity of Method 2 in terms of the m-averaging procednre and the 

effective leakage matrix 

5.1. Sensitivity of Method 2 to details of the m-averaging procednre 

In order to study the sensitivity of Method 2 to the details of the m-averaging procedure, 
as described in Section 12.21 we hrst used the Method 2 code to £t the four different sets 
of m-averaged power spectra iS2010_66a through iS2010_66d (cf. Table [2]). From these 
fits we obtained four tables of htted /- and p-mode parameters: frequencies, linewidths, 
amplitudes, line asymmetries, and their associated uncertainties. From these four tables of 
htted parameters, we collected the frequencies and their uncertainties into four frequency 
tables, J^2010_66a through J^2010_66d, which we then inter-compared on a mode-by-mode 
basis. As one example of these four frequency tables, our table J^2010_66a covered the 
degree range of 0 to 1000, the radial order range of 0 to 29, and the frequency range of 965 
to 7000 /iHz. Within these ranges of /, n, and u we were able to obtain a total of 12,359 sets 
of converged frequencies and frequency uncertainties. 

We note that we have limited our sensitivity study to the investigation of the impact of 
changes in the m-averaged spectra upon the fitted frequencies because the estimates of the 
mode frequencies are an essential part of any structural inversion. 


5.1.1. The influence of the weighting of the un-averaged spectra 

In order to study the possible inhuence upon the fitted frequencies of the weighting 
of the 2/ -|- 1 un-averaged spectra at each degree in the computation of the m-averaged 
power spectra, as described in Section [2T21 we aligned our frequency tables J^2010_66a and 
J^2010_66b on a mode-by-mode basis, and we then subtracted the frequency of each mode in 
J^2010_66a from the corresponding frequency in J^2010_66b. The raw frequency differences 
that resulted, are shown in the two upper panels of Figure HDl where they are plotted, in 
the sense of “weighted” minus “unweighted”, as functions of both frequency and degree. 
Inspection of these two panels illustrates that the largest frequency differences corresponded 
primarily to the high-frequency portions of the ridges for degrees less than 350. 

This concentration of the largest raw frequency differences at the higher frequencies is 
illustrated numerically in section Cl of Table [5l where we show in row 1 that the average 
magnitude of the differences for the frequency range u < 7000 /rHz is six times larger than 
the average magnitude of the cases for which u < 4500 /rHz. On the other hand, when we 
normalized these raw frequency differences by dividing each of them by its formal error, as is 





30 


shown in the lower left panel of Figure [101 we found that only 4.6 % of them were statistically 
signihcant. This is also demonstrated in row 2 of Table |5l where we show that the average 
magnitude of the normalized frequency differences was about 1 a for both frequency ranges 
z/ < 7000 /iHz and z/ < 4500 /rHz, respectively. 

Because most structural inversions have been limited to frequencies that have been less 
than 4500 /iHz, we also wanted to study the radial distribution of the normalized frequency 
differences for the cases that were below this limit. We found that 94.5 % of the cases that 
remained in this restricted frequency range were within the range of ±3 a. Furthermore, when 
we plotted the remaining normalized frequency differences as a function of the fractional 
inner turning-point radius of the corresponding modes, as shown in the lower-right panel 
of Figure [TUI we found that the majority of the cases that were the most signihcant were 
concentrated in the solar convection zone. In fact, very few of the normalized frequency 
differences exceeded ±3 a inward of the base of the convection zone. Overall, Figure [TUI 
indicates that the weighting of the un-averaged power spectra in the computation of the 
m-averaged power spectra had a very modest effect upon the htted frequencies. 


5.1.2. Influence of correction of frequency-splitting coefficients 

In contrast to the situation that we just described in Section [5. 1.11 for the construction 
of the m-averaged power spectra using both unweighted and weighted tesseral and sectoral 
power spectra, the correction of the un-corrected, non-n-averaged frequency splitting co¬ 
efficients for the distortions introduced by the latitudinal differential rotation resulted in 
highly signihcant frequency diherences for those cases that are used in structural inversions, 
as we have demonstrated in Figure [TTJ In the upper-left panel of this hgure, we show the 
frequency dependence of the un-normalized frequency diherences in the sense of “corrected” 
minus “un-corrected” that resulted when we subtracted the frequencies in table J^2010_66c 
from those in table J^2010_66d. In contrast to the upper-left panel of Figure [IHl the ma¬ 
jority of these frequency diherences for the modes having z/ < 4500 /rHz were negative and 
the majority of the high-frequency diherences were positive. In another contrast with the 
upper-left panel of Figure [TOl there is a strong frequency dependence of the high-frequency 
diherences shown in the upper-left panel of Figure [TTl In a third contrast with the situation 
shown in Figure [TUI fhe upper-right panel of Figure [TT] shows that the positive frequency 
diherences were not restricted to the lower-degree modes alone, but instead were present 
for all modes having degrees greater than 200. The negative frequency diherences that are 
shown for 2500 < u < 4500 /rHz in the upper-left panel of Figure [TT] corresponded to the 
rz, = 0 through n = 3 ridges. These negative diherences indicate that the frequencies that 




-31 - 


were computed from the m-averaged power spectra that were computed using the corrected 
splitting coefficients were systematically smaller than were the frequencies for the same ridges 
that were £t to the averaged spectra that were generated using the un-corrected splitting 
coefficients. These systematic frequency differences resulted from the fact that the corrected 
splitting coefficients in this freqnency range were generally smaller than their un-corrected 
counterparts. 

When we normalized these raw frequency differences by dividing each of them by the 
formal uncertainty of the difference, we fonnd that the most signihcant normalized differences 
corresponded to modes which have frequencies between 1800 and 4500 /rHz, as is shown here 
in the lower-left panel of Figure [TT] In particular, the normalized frequency differences were 
highly signihcant for the n = 0 through n = 2 ridges. In the contrast to the situation 
for the raw frequency differences that are shown in the upper-left panel of this figure, the 
normalized frequency differences of the high-frequency modes were mainly less than ±3 a. 
When we restricted these normalized frequency differences to those having frequencies less 
than 4500 /iHz and plotted the remaining cases as a function of the fractional inner tnrning- 
point radius of the modes, as shown in the lower-right panel of Figure [HI we found that 
the cases with the highly-signihcant negative ratios were all concentrated in the shallow 
subphotospheric layers. 

The statistics of both the raw and normalized frequency differences are listed in section 
C2 (i.e., rows 3 and 4) of Table O Here, we show in row 4 that the average magnitude of the 
normalized frequency differences for the frequency range u < 7000 /iHz was equal to 4.9 a, 
while the average magnitude of the normalized freqnency differences for the cases for which 
z/ < 4500 /iHz was equal to 7.8 a. Furthermore, this average of 7.8 a is over 79 times its 
standard error away from zero. Clearly, the correction of the frequency-splitting coefficients 
that are used in the generation of m-averaged power spectra is an essential step. 

A previous comparison of the sets of frequencies computed using both raw and corrected, 
n-averaged splitting coefficients had shown small enough frequency differences that led us to 
believe that the rather time-consuming correction of more than 12,000 raw, non-n-averaged 
splitting coefficients would not be necessary. However, in response to the suggestion of 
the anonymous referee that we investigate the effects of making such corrections, we found 
signihcant diherences in the sets of frequencies compnted using the raw and corrected, non- 
n-averaged splitting coefficients, as is demonstrated here in the fourth row of Table [5l (i.e., 
the second row of Section C2). 
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5.1.3. Influence of non-n-averaged frequency-splitting coefficients 

The influence upon the fitted frequencies of the use of non-n-averaged frequency-splitting 
coefficients in the generation of m-averaged power spectra is illustrated in Figure [121 The 
raw frequency differences, in the sense of “non-n-averaged” minus “n-averaged”, are shown 
as functions of frequency and degree in the upper-left and upper-right panels, respectively. 
These frequency differences resulted when we subtracted the frequencies in table J^2010_66b 
from those in table J^2010_66d. The principal difference between these raw frequency dif¬ 
ferences and those shown in the upper panels of Figure [H] is the absence of the strong, 
wave-like frequency dependence of the high-frequency cases that was seen at the right side 
of the upper-left panel of Figure [TTl It was this wave-like behavior of the high-frequency 
differences that was also visible as the series of shifted peaks that are visible in the upper- 
right panel of Figure [TT] Hence, the absence of such a similar wave-like variation in the 
high-frequency differences that are shown in the upper-left panel of Figure [12] resulted in the 
absence of a similar series of shifted peaks in the upper-right panel of Figure [T^ In turn, the 
absence of this series of positive peaks in the upper-right panel of Figure \n\ meant that the 
majority of the high-degree frequency differences that are shown in the upper-right panel of 
Figure [12] were negative. 

The frequency dependence of the normalized frequency differences that resulted from 
the use of the non-n-averaged frequency-splitting coefficients is shown in the lower-left panel 
of Figure [12] Close comparison of this panel with the corresponding panel of Figure [TT] 
illustrates that both sets of normalized frequency differences had a very similar frequency 
dependence. As was the case in Figure [TTl the most significant of these normalized frequency 
differences were negative and corresponded to z/ < 4500 pHz. 

The dependence of the subset of the normalized frequency differences for which v < 
4500 pHz upon the fractional inner turning-point radius is shown in the lower-right panel 
of Figure [T2j As was the case in Figure [TTl the most significant subset of these normalized 
frequency differences was concentrated in the outer portion of the convection zone. Very few 
of the normalized frequency differences that had inner turning-point radii inward of 0.875 Rq 
were outside the range of ±3 a. 

The statistics of both the raw and normalized frequency differences are given in section 
C3 (i.e., rows 5 and 6) of Table [5] Here, we can see in row 6 that the average magnitude 
of the normalized frequency differences for the frequency range z/ < 7000 /iHz was equal 
to 5.8(7, while the average magnitude of the normalized frequency differences for the cases 
for which z/ < 4500 pHz was equal to 10.2 a. Comparison of both averages with those for 
the same frequency range in row 4 of Table [5] shows that the effects of using non-n-averaged 
splitting coefficients were only slightly more significant than were the effects of the correction 
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of the raw frequency splitting coefficients. Other than this minor difference, the inffuence of 
the non-n-averaged splitting coefficients upon the htted frequencies was very similar to that 
of the correction of the raw frequency splitting coefficients themselves. 

We have just shown in both Figure [12] and in section C3 (i.e., rows 5 and 6) of Table [5] 
that the use of corrected, non-n-averaged splitting rather than the corrected, n-averaged 
splitting coefficients in the generation of the m-averaged power spectra produced signihcant 
frequency differences for the low-order ridges. In order to determine which ridges were 
most sensitive to the changes in the splitting coefficients due to the use of non-n-averaged 
coefficients, we compared the differences between the frequencies in Tables J^2010_66b and 
J^2010_66d with the differences in the two sets of corrected splitting coefficients themselves, 
and we found that the use of the n-averaged splitting coefficients was equivalent to using 
only the splitting coefficients of the higher-order ridges in the generation of the m-averaged 
spectra. The use of the non-n-averaged splitting coefficients for the low-order ridges is what 
resulted in the introduction of the signihcant, systematic frequency differences. 


5.2. Sensitivity of Method 2 to an increase in the widths of leakage matrix 

peaks 

In addition to studying the inhuence on the htted frequencies of the details of the 
computation of the m-averaged power spectra, we also investigated the ehects of an increase 
in the width of the Gaussian approximation that we employ to represent the peaks of the 
leakage matrices (cf. Section 14.ip . To this end, we rehtted the m-averaged spectral set 
i52010_66b with the Method 2 code using a Gaussian approximation the width of which 
has been increased artihcially by 18%. This htting run resulted in the frequency table 
J^2010_66b.lkm. By subtracting the frequencies in table J^2010_66b from those in table 
J^2010_66b.lkm we obtained raw frequency diherences the frequency and degree dependencies 
of which are shown, in the sense of “wider Gaussian” minus “original Gaussian”, in the upper- 
left and upper-right panels, respectively, in Figure HS] Comparison of these two panels with 
the upper panels of Figure [10] shows that these raw frequency differences were similar to those 
shown in Figure [TO] with the main differences being an increase in the number of cases at 
both high frequencies and moderate degrees for which Az/ was positive and a corresponding 
decrease in the number of cases for which Az/ was negative. Further comparison of the 
two upper panels of Figure [T^ with the upper panels of both Figures [H] and [I2] shows that 
the increase in the width of the leakage matrix peak did not result in systematic frequency 
differences for the low-order ridges. 

The normalized frequency differences that resulted from the use of the different approx- 
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imations to the peaks of the leakage matrices are shown in the lower-left panel of Figure [131 
Comparison of this panel with the lower-left panel of Figure [TOl conhrms the absence of any 
signihcant frequency differences in both hgures. The dependence of the subset of the normal¬ 
ized frequency differences for which v < 4500 /iHz upon the fractional inner turning-point 
radius is shown in the lower-right panel of Figure [131 This panel demonstrates that, with 
the exception of a small concentration of signihcant differences just below the photosphere, 
there is no obvious radial dependence to the remaining differences that exceeded ±3 a. 

The statistics of the frequency differences that are shown in Figure [T^ are summarized 
in section C4 of Table [5l Here, row 7 shows that the average magnitude of the raw frequency 
differences was 0.149 /iHz for the cases for which v < 7000 /iHz and dropped to 0.060 /iHz 
for the cases for which v < 4500 /iHz, while row 8 conhrms our statement above that these 
frequency diherences were not signihcant, since the average ratio was equal to 1.00 over the 
entire frequency range, and was equal to only 1.37 for the low-frequency subset. 


6. Sample results from Method 2 and comparison with Method 1 

Before we shall present sample results from Method 2 in Section 16.21 we will demon¬ 
strate in Section 16.11 the substantial improvements that Method 2 makes in the frequencies, 
linewidths, and amplitudes that we generated with this method by comparing them with cor¬ 
responding quantities that we generated using our Method 1. For a description of Method 1 
we refer the reader to Appendix 


6.1. Comparison of results from Methods 1 and 2 

One key diherence between Method 1 and our Method 2 is the fact that Method 1 is 
restricted to use only symmetric Lorentzian prohles, while Method 2 includes the option 
of using either symmetric or asymmetric prohles. A second key difference between the 
two methods is the use of both narrow and wide htting ranges in Method 1, while such 
wide htting ranges are never employed in Method 2. Instead, in Method 2 we employ the 
procedure described in Section 14.31 to determine the width of each htting box individually. 
As a result, the widths of the htting boxes vary systematically as functions of both the degree 
and the frequency in Method 2. 

In the three left-hand panels of Figure [TH, we compare the htted prohles that Method 1 
produced with segments of three diherent spectra from our m-averaged spectral set iS2010_66a. 
In the three right-hand panels of the same hgure we compare the prohles that Method 2 pro- 
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duced for the same three spectral segments using the asymmetric prohle of iNigam fc Kosovichev 
fll998l) . The normalized comparison of Method 1 and Method 2 £t results shown in Table [6] 
for the same three modes whose fits are shown in Figure [TTI demonstrates that the Method 2 
frequencies and linewidths do agree more closely with the Method 1 results that were com¬ 
puted using the narrow htting range than they do with the Method 1 results that were 
computed using the wide htting range. Table |6] also shows that, with the exception of the 
frequency that was computed for the {n,l) = (2, 70) mode using the narrow htting range in 
Method 1, all of the other frequencies and linewidths are systematically diherent between 
Method 1 and Method 2. 


In Figure [15] we illustrate the improvements that Method 2 made in both the formal 
frequency and linewidth uncertainties relative to the corresponding quantities that were 
computed using Method 1. In the upper-left panel of Figure [T5] we show the frequency 
dependence of the ratios of the Method 1 and Method 2 frequency uncertainties for the cases 
in which the narrow htting range was employed in Method 1, while we show the frequency 
dependence of the similar ratios for the cases in which the wide htting range was employed 
in the upper-right panel. The average ratio of the two sets of uncertainties in the upper-left 
panel of this hgure was 36.1 ± 1.5, while the average ratio of the two sets of uncertainties in 
the upper-right panel was 9.7 ± 0.1. 

The frequency dependence of the ratios in the Method 1 and Method 2 linewidth un¬ 
certainties are shown in the lower-left panel of Figure [T^ for the cases in which the narrow 
htting range was used in Method 1, and we show the frequency dependence of the ratios of 
the Method 1 and Method 2 linewidth uncertainties for the cases in which the wide htting 
range was employed in the lower-right panel of the same hgure. The average ratio of the 
two sets of linewidth uncertainties in the lower-left panel was 96.2 ± 4.3, while the average 
ratio of the two sets of linewidth uncertainties in the lower-right panel was 12.6 ± 0.1. Com¬ 
bining both the narrow and wide htting range ratios in overall sets of uncertainty ratios, we 
found that Method 2 reduced the frequency uncertainties by an average ratio of 15.2 ± 0.3 
and it reduced the linewidth uncertainties by an average ratio of 29.8 ± 0.9. These sizeable 
reductions in both the frequency and linewidth uncertainties testify to the importance of 
including more than a single peak in the theoretical prohle that is used to ht both modes 
and ridges of power. 

In the upper two panels of Figure [16] we show the improvements that Method 2 in¬ 
troduced into the frequencies of the n=0 ridge in comparison with the frequencies com¬ 
puted using Method 1. We demonstrate these improvements by comparing both sets of 
htted frequencies with theoretical frequencies that have been compute d from Model S o f 
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Specifically, in both of these panels we show the freqnency dependence of the freqnency dif¬ 
ferences Au = where refers to either our Method 1 or Method 2 htted fre¬ 

quencies and where is the corresponding theoretical frequency computed using Model S. 
The Au values for the Method 1 hts are shown as the black diamonds, while the Au values 
for Method 2 are shown as the red diamonds. The upper-left panel is for the portion of the 
n = 0 ridge that was originally £t using the narrow fitting range in Method 1, while the 
upper-right panel is for the frequency range of the same ridge that was initially £t using the 
wide htting range in Method 1. The upper-left panel shows that the Method 2 frequencies 
are much smoother than the Method 1 frequencies and that they also agree more closely 
with the theoretical frequencies. The upper-right panel shows that, with the exception of a 
range of frequencies around 2700 pHz, the Method 2 frequencies agreed more closely with 
the theoretical frequencies than did the Method 1 frequencies. 

In the lower-left panel of Figure [16] we show the improvements that Method 2 introduced 
into the linewidths of the n = 0 ridge. Since there are no theoretical linewidths that we can 
compare the fitted linewidths with, we have simply plotted both sets of fitted linewidths 
as a function of frequency in this panel. This comparison demonstrates that the Method 2 
linewidths (shown as the red diamonds) varied in a much smoother manner with increasing 
frequency than did the Method 1 linewidths (shown as the black diamonds). It also demon¬ 
strates that for the majority of the modes shown, the Method 2 linewidths were smaller than 
were the corresponding linewidths computed using Method 1. 

In the lower-right panel of Figure [12] we extend our comparison of the Method 1 and 
Method 2 linewidths to the entire set of 30 ridges that were ht using Method 2 and using both 
fitting ranges in Method 1. Here we have plotted the frequency dependence of the ratios of 
the Method 2 linewidths divided by the Method 1 linewidths. This panel clearly shows that 
for all of the modes with frequencies below the acoustic cut-off frequency the fits computed 
using Method 2 produced linewidths that were smaller than the corresponding Method 1 
linewidths by factors ranging between 2 and 5. The smaller linewidths all correspond to 
longer estimated lifetimes for these modes. 

in order to show the entire frequency range spanned by both htting ranges with a 
common scale for the y-axis” 

In the upper-left panel of Figure [T7] we have combined both portions of our comparisons 
of Method 1 and Method 2 frequencies with the Model S frequencies for the n = 0 ridge into 
a single panel in order to show the entire frequency range spanned by both fitting ranges 
with a common scale for the vertical axis to illustrate yet another important improvement 
of Method 2 over Method 1 - namely the complete absence in the Method 2 frequencies 
of the substantial discontinuity that is present in the Method 1 frequencies at precisely the 
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frequency where the transition was made from the use of the narrow htting range to the 
use of the wide htting range in that method. In the upper-right panel of Figure [T7] we 
show a similar comparison of Method 1 and Method 2 frequencies with the corresponding 
Model S frequencies for the n = 1 ridge. As with the upper-left panel, there is a pronounced 
discontinuity in the Method 1 frequencies at the location of the switch-over from the narrow 
htting range to the wide htting range that is absent in the Method 2 frequencies. 

In addition to the introduction of jumps into the htted frequencies, the switch-over in 
the width of the htting ranges that were employed in Method 1 also caused a corresponding 
discontinuity in the Method 1 linewidths, as we show here in the two middle panels of 
Figure fT71 The jumps in the linewidths that are shown in black in these two panels occurred 
for exactly the same modes for which the frequencies also exhibited the jumps in the two 
upper panels. In contrast to the Method 1 linewidths, the Method 2 linewidths exhibited 
no such discontinuities, as we show in red in the middle panels of Figure [T71 As we also 
showed in the lower panels of Figure fT6l the Method 2 linewidths agreed more closely with 
the Method 1 linewidths that were computed using the narrow htting range than they did 
with the Method 1 linewidths that were computed using the wide htting range. 

We also demonstrate in the two lower panels of Figure [T7] that the use of the single¬ 
peak prohle in Method 1 introduced discontinuities in the amplitudes (as shown in black), 
or power densities, that disappeared when the multiple-peak prohle of Method 2 was used 
instead (as shown in red). Not only did the single-peak prohle introduce discontinuities into 
the amplitudes, it also caused the amplitudes that were computed using the narrow htting 
range to show more scatter as a function of degree along the ridges to the left of the jumps 
than was the case with our Method 2 amplitudes. 

In Figure [TS] we demonstrate that the replacement of the narrow htting range with the 
wide htting range in Method 1 did not just introduce distinct discontinuities in the htted 
frequencies, linewidths, and amplitudes for the n = 0 and n = 1 ridges, but rather introduced 
similar discontinuities in all three quantities for nearly all of the ridges that we ht. We found 
it useful to measure these discontinuities in terms of the normalized discontinuities dehned 
as 

(40) 

where T can be either frequency z/, linewidth ta, or else the logarithm of the amplitude A, 
and where 

6T{h, k) = AT{k) - AT{h) = r fi*(/2) - , (41) 

with li being the degree of the highest-degree mode that we ht using the narrow htting range 
for a given ridge, with I 2 being the degree of the lowest-degree mode that we ht using the 
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wide fitting range for the same ridge, and with T®* and denoting, respectively, the 

htted and the corresponding seed value of T, and where is the point-to-point scatter in 
the AT values for the same ridge, as computed from equation fIBip using 

V = AT/Al, (42) 

with AT = T®* — The use of the differences T®* — T®®®'^ rather than the htted values, 

T^*, themselves in the above equations 0411) and 042p has the advantage that the changes 
AT /Al are approximately zero along the ridge under consideration and, thus, gives rise to 
an unbiased value of the scatter. Sat, as pointed out in Appendix IBll We note that the 
values of T®®®*^ that we have employed are smoothly varying with degree, and we also note 
that we used a different pair of li and I 2 values for each of the analyzed ridges, depending 
on the degree at which the narrow htting range is replaced with the wide htting range in 
Method 1. 

In the upper-left panel of Figure fT8l we show the normalized discontinuities, (5z/(/i, I 2 ) / Sau, 
in both the Method 1 and Method 2 frequencies in dependence of the radial order n. These 
normalized discontinuities were computed from equations fl40p through 0421) using T = z/, 
with u denoting the htted frequency from either Method 1 or Method 2. The hereto required 
values of the seed frequency, were taken from a suitable seed table (cf. Section ITTSp . 

In order to help indicate which of the normalized frequency discontinuities represent sta¬ 
tistically signihcant values, we have included the two dashed green lines at ratios of ±3. 
Inspection of the upper-left panel of Figure [TS] shows that for Method 1 (black diamonds) 
14 of the 30 normalized frequency discontinuities exceeded ±3 a. By contrast, of the 30 nor¬ 
malized frequency discontinuities that we generated using Method 2, and which are shown 
as the set of red triangles in the upper-left panel of Figure [T51 none of them exceeded ±3 a. 
Clearly, the multiple-peak prohle that we are using in Method 2 has allowed us to overcome 
the discontinuities that we were forced to introduce into the htted frequencies when we were 
using the single-peak prohle of Method 1. 

The introduction of the jumps into the Method 1 frequencies that corresponded to 
the change-over in the htting range that we had to employ with Method 1 meant that those 
frequencies could not be employed in structural inversions since these jumps would introduce 
unphysical oscillations in the radial prohle of any quantity determined in such inversions, e.g., 
internal sound speed. We have re-plotted the normalized Method 1 frequency discontinuities 
that are shown in the upper-left panel of Figure [TH] as a function of the inner turning-point 
radius of the li mode in the upper-right panel of Figure 0 It is clear from this panel that 
one group of signihcant Method 1 frequency discontinuities corresponds to modes that span 
the lower portion of the convection zone, while the other group of statistically signihcant 
discontinuities corresponds to modes that span the sub-surface shear layer. On the other 
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hand, the fact that none of the 30 normalized frequency discontinuities computed using 
Method 2 exceeded ±3 a means that this method is capable of producing frequencies that 
can be employed in successful structural inversions, as we will show later in section [3 

For the linewidths, as shown in the lower-left panel of Figure flSl 23 of the 30 normalized 
Method 1 discontinuities, 6 w{li, l 2 )/^Aw, exceeded ±3 a, while only the n = 0 value exceeded 
±3 a for Method 2. For the logarithm of the amplitudes, as shown in the lower-right panel of 
Figure [T51 21 of the 30 normalized Method 1 discontinuities, 6 \ogA{li,l 2 )/^AiogA, exceeded 
±3 a. However, none of the normalized Method 2 discontinuities exceeded ±3 a. Clearly, 
Method 2 is able to produce both linewidths and amplitudes which do not exhibit the 
systematic discontinuities that Method 1 introduced into both quantities. 

In addition to all of the discontinuities in the frequencies, linewidths, and amplitudes. 
Method 1 also introduced similar discontinuities into the frequency uncertainties, as we 
show in both panels of Figure [191 In these panels we show the frequency dependence of the 
logarithm of the Method 1 frequency uncertainties that were computed from the hts to the 
m-averaged spectral set iS2010_66a as the black curves. The left-hand panel is for the n = 1 
ridge, while the right-hand panel is for the n = 2 ridge. There are clearly discontinuities in 
both curves at the locations of the vertical dashed lines, which are located at the frequencies 
of the transition from the narrow fitting range to the wide fitting range for both ridges. By 
contrast, the corresponding frequency uncertainties that we computed for the same power 
spectra for both of these ridges using Method 2 are shown as the two red curves. It is clear 
that neither of these curves contains a discontinuity at the location of the vertical lines. 


6.2. Sample results from Method 2 

The results that we will now describe were obtained from fitting the m-averaged spec¬ 
tral sets iS2010_66a and iS2010_03 using Method 2. These two spectral sets have been com¬ 
puted from the 66-day long observing run 7^2010_66 and from the 3-day long observing run 
7?.2010_03, respectively. The rationale behind the htting of spectra from both a 66-day long 
time series and a 3-day long time series is to investigate the impact of spectral resolution 
upon the performance of Method 2. In each of these two sets of m-averaged spectra we 
considered modes in the range 0 < / < 1000, 0 < n < 29, 900 < v < 7000 /iHz. In this 
range of frequencies, radial orders, and degrees the seed table that we had available to us 
contained a total of 12,533 entries. 

In the low-degree, low-frequency portion of the dispersion plane, the mode amplitude 
becomes comparable to the background noise, and at the same time the mode linewidth 
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becomes smaller than the spectral resolution. This can easily lead to the confusion of a 
noise spike with a modal peak. By a careful inspection of the quality of the fits we therefore 
set, for each radial order, n, a lower degree limit, below which we did not try to fit a 
mode for that ridge. In a similar analysis of the high-degree, high-frequency portion of the 
dispersion plane we defined, for each ridge, an upper degree limit, /mix, above which we did 
not try to fit a mode for that ridge. Based upon this approach we came up with a total of 
12,359 modes to be £t for the 66-day spectra, and a slightly smaller total of 12,242 modes to 
be £t for the 3-day spectra because of the lower signal-to-noise ratios present in this set of 
spectra. We found that for all of the modes that we attempted to fit in the case of the 66-day 
spectra and the 3-day spectra, respectively, the Method 2 code converged on a solution. As 
a result, we ended up with a total of 12,359 sets of fitted mode parameters for the 66-day 
spectra, and a total of 12,242 sets of htted mode parameters for the 3-day spectra. This 
impressively demonstrates that Method 2 is working in a stable manner over wide ranges of 
frequency, degree, and spectral resolution. 


We collected the 12,359 frequencies and their corresponding uncertainties that we gener¬ 
ated with Method 2 from the iS2010_66a m-averaged power spectra into the frequency table 
J^2010_66a (cf. Section [5Tll . The dispersion plane coverage of this set of frequencies is illus¬ 
trated in Figure [201 For the reasons explained in Section 14.61 we employed the symmetric 
Lorentzian prohle for fitting the pseudo-modes. For all of the cases for which the Lorentzian 
prohle was used, the line asymmetry was defined to be zero, i.e., B = 0. Those cases are 
shown as the red diamonds in Figure l20l For all of the other cau s es, wh ich are shown as the 
black diamonds, the asymmetric profile of lNigam fc KosovichevI fll998l) was employed. 


In addition to the frequencies that we computed from the iS2010_66a m-averaged power 
spectra using Method 2, we show the corresponding linewidths and line asymmetries in 
the upper-left and lower-left panels of Figure EH respectively. The upper-left panel shows 
that the linewidths of the higher-order ridges exhibited either local or global maxima for 
frequencies between 5800 and 6000 /iHz before decreasing slightly at higher frequencies. The 
frequencies where the linewidths of the ridges exhibited either local or global maxima are 
shown as the green open circles in Figure l20] that are connected by the solid green line. All of 
these frequencies can be seen to lie well above the range of frequencies where the fitting prohle 
was switched from an asymmetric to a symmetric prohle, which was 4638 < u < 4842 /iHz. 


The maxima in the linewidt hs are very similar t o peaks in the high-frequency linewidths 
that were shown in Figure 4 of iDuvall et ah! fll99ll) . If we compare the frequencies of the 
peaks in the linewidths in the degree range that they employed (i.e., 11 < / < 150) with the 
frequencies as given by the green line in Fig ure l20l we see that our frequencies are about 
200 /iHz larger than the frequencies given bv iDuvall et ah! fjl99ll ). This small diherence is 














41 


most likely due to the different levels of solar activity in the two observing epochs. The 


Duvall et al.l fllQQlh observations were acquired in late 1987, while the 66-day data set that 


we have employed, viz. iS2010_66a, was acquired in mid-2010. We hnd it intriguing that the 
solid green line that connects the maxima in the linewidths in Figure 120) loo ks very similar to 
the de gree dependence of the chromospheric modes as shown in Figure 5 of lUlrich fc Rhodes 
(119771) : however, the frequencies of these linewidt h maxima are all much closer to the fre¬ 
quency of the second chromospheric mode, which lUlrich fc RhodesI (119771 ) computed to be 
5570 /iHz than they are to the frequency of the hrst mode, which they computed to be 
4138 /iHz. Furthermore, there is no evidence in Figure [21] for the hrst chromospheric mode 
in the linewidths, nor is there any evidence in Figure [20] for either chromospheric mode in 
the frequencies themselves. Since there is no other evidence for either chromospheric mode 
in any of our frequency tables, the origin of these local maxima in the linewidths is not clear, 
but is probably related to the transition from standing waves below the acoustic cut-off fre¬ 
quency to pseudo-modes above. The lower-left panel of Figure [21] shows that most of the line 
asymmetries were negative, while they were zero by construction for those cases for which 
the symmetric Lorentzian prohle was used in Method 2. 


The amplitudes, or power densities, as is shown in the upper-right panel of Figure EH 
peak very close to 3000 /rHz, in agreement with all previous studies of the solar oscillation 
amplitudes. In contrast, the power values that are shown in the lower-right panel go through 
their peak at about 3300 /rHz because the power values are computed from the products 
of the amplitudes and the linewidths and the increase in the linewidths with increasing 
frequency partially compensates for the decreasing amplitudes. As we summarized at the 
beginning of this section, we also obtained similar sets of linewidths, line asymmetries, 
amplitudes, and power values for the set of 3-day power spectra in addition to the results 
we have shown in this section from the 66-day power spectra. 


The smoothness in both the htted frequencies and linewidths that resulted from the 
hts to the iS2010_66a m-averaged power spectra using Method 2 is illustrated in Figure [22] 
using the example of the n = 14 ridge. The upper-left panel shows the degree dependence 
of the frequencies for that ridge, and it also includes the smooth £t to those frequencies as 
the red curve. The lower-left panel shows the degree dependence of the frequency slope, 
Au/Al, that was computed as described in the hgure caption, and it also includes, as the 
blue curve, the derivative with respect to degree of the red curve shown in the upper-left 
panel. The fact that the derivative with respect to degree (blue line in the lower-left panel) 
of the smooth curve (red line in the upper-left panel) htted to the frequencies is an excellent 
representation of the numerical values of Au/Al (black diamonds in lower-left panel) at all 
of the degrees shown, impressively demonstrates the smoothness of the htted frequencies 
(black diamonds in upper-left panel). The upper-right panel of Figure [22] shows the degree 
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dependence of the linewidths for the n = 14 ridge, and it also includes the smooth £t to 
those linewidths as the red curve. This panel shows that this ridge was one of the ridges for 
which the linewidth exhibited a local maximum before decreasing toward the high-degree 
ends of those ridges as we just described in Figure [211 The lower-right panel of Figure 1221 
shows the linewidth slopes, Aw/Al, which were computed as described in the hgure caption, 
as a function of degree, and it also includes, as the blue curve, the derivative of the smooth 
curve with respect to degree that is shown as the red curve in the upper-right panel. As was 
the case for the frequencies, the linewidths are also smooth enough that the derivative of the 
smooth fitting curve provides an excellent representation to the values of Aw/Al. 

As we noted in Section 14.21 the derivatives with respect to degree of the frequency, 
the linewidth, and the amplitude are an essential ingredient of Method 2. We compute 
these derivatives on a ridge-by-ridge basis by fitting, as a function of degree, higher-order 
Chebyshev polynomials to the frequencies, linewidths, and amplitudes, respectively, as they 
have been determined with Method 2, and by then differentiating the resulting polynomials 
with respect to degree. The red curves in the upper panels of Figure [22] are examples of 
such hts to the frequency and linewidth, respectively, while the blue curves shown in the 
lower panels are examples of the hrst derivative of the htted polynomials. Because higher- 
order polynomials are involved, this approach is prone to deliver wildly oscillating htting 
curves. Therefore, the degree of the polynomials has to be selected rather carefully, and the 
least-squares hts of the polynomials need to be constrained as to avoid oscillating curves. 


Helioseismic inversion for solar structure 


In this Section we will present the results of an inversion for the radial structure of 
the Sun using a subset of the frequencies that are contained in our table J^2010_66a. We 
have already mentioned in Section 16.21 that this entire table contains a total of 12,359 sets 
of htted modal parameters, and we showed the dispersion plane coverage of this complete 
set of frequencies in Figure [201 This is the largest mode-set to have ever been ht thus 
far in helioseismology. To carry out our structural inversion, we selected a subset of 6,366 
frequencies which contained degrees that ranged from 0 to 1000, radial orders from 0 to 
29, and frequencies from 969 to 4500 /iHz. Since this frequency upper limit of 4500 pHz is 
less than the frequency of the lowest symmetrical ht, which was 4638 /iHz for the n = 17 
ridge, all of the frequencies in this subset were computed using the asymmetric prohle of 
Nigam fc KosovichevI fjl998[) . By using the Optimally Localized Averaging (OLA) technique 
we inverted this subset of modes to determine the spherically symmetric structure of the Sun. 
We limited the frequency range of the subset of modes by 4500 /iHz because this threshold 
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is well below the acoustic cut-off frequency which separates the p-modes from the pseudo¬ 
modes and which is about 5000 /iHz. We note that the choice of the upper frequency limit 
of the mode set to be inverted is not critical since modes with frequency above 4000 /iHz 
do not contribute signihcantly to the inversion results due to their large errors. Details of 
the inversion procedure, i nclud ing calculation of the sensitivity kernels and test results, are 
presented by iKosovichevI (119991) . 


While we had to select from our full mode-set a suitable subset of frequencies for our 
structural inversion, we note, however, that tables of htted mode parameters covering a simi¬ 
lar range of degrees and frequencies as our full mode-set presented here, are an indispensable 
prerequisite for the study of temporal changes in the sensitiviti es of the mode parameters to 


corresponding changes in the levels of s olar activity (see, e.g.. iRabello-Soares et ahll2008b 


Rhodes et al. 

2011; 

Rabello-Soares 

2011) 


7.1. Analysis of systematic errors in fitted frequencies and associated 

uncertainties 

In an inversion for solar internal structure both the frequencies and the associated 
uncertainties of the solar oscillation modes are employed as input parameters. Therefore, 
systematic errors in these input parameters may affect the resulting inferences of the Sun’s 
internal structure. As we have explained in Section 13.2.31 such systematic errors may arise 
from instrumental effects. However, they also are likely to arise from the methodology em¬ 
ployed for htting the spectra. For example, the use of an inadequate htting box around the 
target peaks, an inadequate treatment of the background power within the selected htting 
boxes, and/or contributions from n- and /-leaks that are not adequately accounted for in the 
htting model prohle, as given in equations fl24|) through fl2^ . may lead to systematic errors 
in both the htted frequencies and the uncertainties thereof. On the other hand, the presence 
of a group of outliers in a table of frequencies and associated uncertainties that is designed 
to be used as input to a solar structural inversion, may result in the appearance of unphys¬ 
ical oscillations in the radial prohles of the resulting sound speed or other thermodynamic 
quantities. Such a set of outliers may also have the unwanted side ehect of preventing a 
stable, regularized solution from being found outside of a very narrow range of the regular¬ 
ization parameter that is being employed in the inversion procedure. In order to avoid the 
appearance of such unphysical oscillations in the inverted sound speed prohle, we developed 
a hve-step procedure that allows us to substantially alleviate the issue of systematic errors in 
both the htted frequencies and their associated uncertainties in tables that will be employed 
as input data sets for structural inversions. The rather high complexity of this procedure 

















44 


mainly results from our desire to keep the approach as objective as possible, while keeping 
unavoidable subjective elements at a minimum. 


The basic idea behind our correction scheme is to compare the observed frequency for 
the (n, /) mode, where the term “observed frequency” refers to our original set of htted 
frequencies, with the theoretical model frequency for the same mode, on a ridge- 

by-ridge basis. For this comparison we employ two sets of frequency differences between 
the observed frequencies, and the model frequencies, viz. the set of unsealed 

differences. 
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and the set of scaled differences 
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In equations fl4^ and (l44)) v™, is the theoretical frequency of mode (n, /) that has been 
computed from Model S of IChristensen-Dalsgaard et al.l (119961) employing the approach de¬ 
scribed bv iKosovichevI (119991) . and is the inertia of mode (n, /), normalized to the inertia 
of mode (u,0). We note that both and qn^i are computed using the same solar model 
that is also used to compute the kernel functions employed in the structural inversion. The 
rationale behind the use of two different sets of frequency differences, Avn.i and di'n.h respec¬ 
tively, is to exploit the fact that, for a given ridge, the run of Az/„ ^ with respect to degree 
differs markedly from that of 5i>n,i- This allows the identiheation of outliers that would 
otherwise remain undetected if only one set of frequency differences is used. The rationale 
behind the dehnition of the scaled frequency differences, 5vn,h is to remove the frequency 
differences between the solar model and the Sun th at come from the dominant near-surface 
effects le.g.. lChristensen-Dalsgaard fc Goughlll984i) by approximating these differences using 
a frequency power-law function scaled with the mode inertia. 


Before we delve into the details of our correction scheme, we briefly summarize its ma¬ 
jor steps. First, for each ridge we ht, as a function of degree I, a Chebyshev polynomial, 
of degree cu to the unsealed frequency differences, Ai>n,h and a Chebyshev polyno¬ 
mial, of degree vr to the scaled frequency differences, Although both and 

U-jr are Chebyshev polynomials of the same kind but generally of different degree, we have 
introduced different symbols for them in order to highlight that they are used for the htting 
of different sets of frequency differences. 

Next, from the two ht polynomials, S^\l) and respectively, we generate, for 

each ridge, two sets of emulated (imitated) frequencies, and respectively. Here, 
corresponds to an observed frequency that would fall exactly onto the htted curve 
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while corresponds to an observed frequency that would fall exactly onto the 
htted curve n^\l). At this point we have hve sets of frequency differences at our disposal, 
namely, 
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from which we identify outlying modes by applying suitable statistical criteria. 

Next, we reht, by means of the WMLTP method employing updated seed tables, all of 
the modes that have been flagged as outliers. 

Finally, we use the rehtted frequencies and associated uncertainties to update all of the 
outlying modes. In some cases we only update the frequencies of those outliers and in other 
cases we update only the associated uncertainties, and in still other cases we update both 
the frequencies and their uncertainties at the same time. We also encountered a few cases, 
however, for which neither the frequency nor the uncertainty needed to be altered. 


7.1.1. Step 1 

In the hrst step of our outlier correction scheme we £t, in the least-squares sense, the 
Chebyshev polynomial, of degree vu to the unsealed frequency differences, Aun,i, as 

dehned in equation (03]), along a ridge of given radial order n, i.e., as a function of degree 1. 
Typical examples are shown here in Figure [23] for the n = 1 and n = 11 ridge, respectively. 
In both panels of Figure [23]the unsealed frequency differences, Ah'n,h are shown as the black 
diamonds, while the htted Chebyshev polynomials, are shown as the red lines. The 

degree w of the polynomials S^\l) depends on the radial order, n, and is determined by trial 
and error. In doing so, two requirements must be considered. On the one hand, the degree 
w needs to be large enough so that the htted polynomial E^\l) provides a reasonably good 
ht to the data. On the other hand, it must be small enough to avoid unrealistic oscillations 
in the htted curve. For the ridges in the range 0 < n < 29 we found values of vj in the range 
1 < -07 < 26 to be useful. In particular, we found that w decreased as the order of the ridge 
increased. Unfortunately, so far we have not succeeded in finding any pre-set criteria which 
would allow the determination of the degree w in an automated fashion. 
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If for any mode along the ridge, n, the criterion 


\Aunj - > 3 \/(A<f)"+ (aS<,">(/))", (50) 

is fulhlled, the mode (n, 1) is flagged as an outlier. Here, Aun^i is the unsealed frequency 
difference, is the uncertainty of the observed frequency S^\l) is the htted 

polynomial, and AS^\l) is the uncertainty thereof. While Au^f is calculated by default 
within the WMLTP htting code as described in Appendix AS^\l) is computed from 

AA£")(/) = (51) 


where is the 100 (1 — 7 / 2 ) percentage point of Student’s t-distribution with 

\ — vj — 1 degrees of freedom, A being the number of data points involved in the £t of 
to the frequency differences, Avn,i- As mentioned in Appendixfor la confidence 
intervals 7 = 0.31731. If we denote the coefficients of the Chebyshev polynomial S^\l) by 
Xi, i = 0,... the variance of S^\l) is given by 


var. 
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EE 
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(52) 


f|NAG Fortran Library. Mark 23ll201ll) . Here, S is the sum of squares, and denotes the 
elements of the inverse of the Hessian Matrix. 


Once an optimum value of oj has been found for a given ridge, the outliers are flagged 
by means of the criterion (l50|) . and a new polynomial S^\l) is ht to the remaining points 
using the same value of uj. By means of this rehtted polynomial an updated set of outliers is 
now identihed, and the entire process is repeated until all of the outliers have been identihed 
for that ridge. 

For further analyses we need to generate a set of emulated frequencies 
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eml 

n,l 
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(53) 


which correspond to observed frequencies that would fall exactly onto the curve for a 

given ridge. Because the theoretical model frequencies, are error-free, the uncertainty 

of is given by 

Au:f = AS^\l), 

with AA^\l) being given in equation flSTl) . 


(54) 
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7.1.2. Step 2 


In the second step of our outlier correction scheme we ht the Chebyshev polynomial, 
n^\l), of degree vr to the scaled frequency differences, 5'Un,h dehned in equation fl44|) . 
along a ridge of given radial order, n, i.e., as a function of degree 1. We note that the outliers 
that were identihed during step 1 are not used in the htting of the n^\l) polynomials 
to the scaled frequency differences, 5i'n,i- Typical examples of hts to the scaled frequency 
differences are shown here in the left-hand panels of Figure [241 for the n = 1 and n = 11 ridge, 
respectively. In both of these left-hand panels of Figure l24l the scaled frequency differences, 
5'UrL,h shown as the black diamonds, and the htted Chebyshev polynomials, n^\l), are 
shown as the red lines. The degree tt of the polynomials n^\l) depends on the radial order, 
n, and is determined by trial and error. For the ridges in the range 0 < n < 29 we found 
values of tt in the range 1 < tt < 29 to be useful. Rather than simply decreasing with 
increasing n, as was the case for vj, vr increased from 19 at n = 0 to 29 at n = 4 before 
decreasing to 1 at n = 29. As with the degree vj, we have not yet been able to hnd any 
approach that would allow us to determine the degree tt of the polynomials n'^\l) in an 
automated manner. 

Next in this second step, we determine, for each ridge of radial order, n, the htting 
residuals, of the ht of n^\l) to the scaled frequency diherences, 6i'ri,h i-e., 

= 6vn,i-n^:^\l). (55) 


We hag those modes (n, 1) as outliers for which the criterion 


> 3 std(r|''^^) 


(56) 


is fulhlled. Here, std(rj"^) denotes the standard deviation of the htting residuals, for 
the ridge of radial order, n. For each ridge, n, the ht of the polynomial, to the scaled 

frequency diherences, is repeated until all outliers have been identihed by means of the 
criterion (1561) . 


For further analyses we will need to generate yet another set of emulated frequencies. 


viz. 
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(57) 


which correspond to observed frequencies that would fall exactly onto the curve that 

has been htted to the scaled frequency diherences, Siy'n,i, for the ridge of radial order, n. 
Because the theoretical model frequencies, are error-free, the uncertainty of is 
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given by 
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Qn,l 



3.5 
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(58) 


with Ani'^\l) being compnted from eqnation (l5T]) with ^m\l) being replaced with n^\l). 


In the ht of the Chebyshev polynomial Sm\l) to the unsealed frequency differences, 
Aun^i, along a given ridge (cf. Section 17. 1.11) as well as in the ht of the Chebyshev polynomial 
n^\l) to the scaled frequency differences, along a given ridge, outliers can only be 

detected in a reliable manner if the rough trendline of the htted curve is known as a function of 
degree a priori. While this basic requirement is fulhlled for the unsealed frequency differences, 
Avn,h because their trendline along a given ridge is only a slowly varying function of degree 
(cf. Figure [23D. it is not fulhlled for the scaled frequency diherences, 5i'n,h as is demonstrated 
here in the left-hand panels of Figure [24l for both the n = 1 and n = 11 ridge. For the 
n = 1 ridge a sharp variation of the trendline with degree is prominent at around I ~ 60, 
as is shown in the upper-left panel of Figure [MJ while for the n = 11 ridge similar sharp 
variations do occur at about I ~ 1 and I ~ 7, respectively, as is shown in the lower-left panel 
of Figure [241 Such sharp variations of the trendline with degree raise the question whether 
or not they are real features. However, this is important to know because suspicious trends 
may reflect some essential variations in the solar structure that are not yet known. 


In dealing with this issue we found it useful to introduce the emulated scaled frequency 
differences, given by 
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which differ from the original scaled frequency differences, merely by the replacement of 
the observed frequencies, in equation fl44]) with the emulated frequencies, , defined 
in equation (l53l) . While for some ridges the degree dependence of the trendline of the curve 
htted to the original scaled frequencies, 5i>n,ii cannot be determined in a unique manner 
mainly because of lurking outliers in the observed frequencies, , the degree dependence of 
the trendline of the emulated scaled frequency differences, given uniquely. Primarily, 

this is because the march of the observed frequencies, , is only a slowly varying function 
of degree without any sharp variations so that the march of the emulated frequencies, 
with degree is given unambiguously for each ridge by means of a straightforward ht of the 
Chebyshev polynomial to the unsealed frequency diherences, Ai'n,i (cf- Section [7.1.ip . 

On these grounds, the degree dependence of the trendline of the emulated scaled fre¬ 
quency diherences, can be considered, for a given ridge, as a reliable template for 








49 


the trendline of the curve htted, as a function of degree, to the original scaled frequency 
differences, 6un^i. Therefore, by a careful identihcation of lurking outliers in the observed 
frequencies, as well as a careful selection of the curve (i.e., Chebyshev polynomial n^\l) 
of degree tt) that is htted to the original scaled frequency differences, along a given 

ridge, the degree dependence of the trendline of this htted curve can be made similar to that 
of the emulated scaled frequency diherences, This is demonstrated here in Figure [24] in 

which we show the degree dependence of the emulated scaled frequency diherences, 
the n = 1 ridge in the upper-right panel, and for the n = 11 ridge in the lower-right panel, 
respectively. A close comparison of the left-hand panels and the corresponding right-hand 
panels of Figure [2l| clearly indicates that for both the n = 1 and n = 11 ridge the trend¬ 
lines of the curves htted to the original scaled frequency diherences, 5vri,i-i as a function of 
degree are in excellent agreement with the corresponding trendlines of the emulated scaled 
frequency diherences, We therefore can safely presume that the sharp variations with 

degree to be seen in the degree dependence of the trendlines of the original scaled frequency 
diherences, 5i'n,h of both the n = 1 and n = 11 ridge are real features. We note that the 
same analysis is also useful for ridges other than the n = 1 and n = 11 ridge. 

At the end of the second step of our outlier correction scheme we have the following three 
diverse sets of frequencies and associated uncertainties at our disposal: the set of observed 
frequencies and the uncertainties thereof, as well as the two sets of emulated 

frequencies and respectively, the uncertainties of which are given by and 

respectively. 


7.1.3. Step 3 

In the third step of our outlier correction scheme we compute from these three sets 
of frequencies and associated uncertainties a total of three sets of normalized frequency 
differences, viz. 
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We flag those modes (n, 1) as outliers for which at least one of the following criteria is fulhlled: 

IWI>3. |Ahj|>3. (63) 

7.1.4- Step 4 

The fourth step of our outlier correction scheme actually consists of two sub-steps. As 
we already have pointed out in Section ing custom-made seed tables are required as input to 
the WMLTP method. This is because the fitting profile, as given by equations fIMl) through 
fl29|) . depends on the initially not very well known variation of frequency, amplitude, and 
linewidth with degree, /. Therefore, in the first snb-step we update those seed tables by 
performing an additional iteration of our fixed-point iteration, as described in Section 14.21 
In the second snb-step we then reht, by means of the WMLTP method, all of the modes that 
have been flagged as ontliers in steps 1 through 3 of our outlier correction scheme, thereby 
employing the updated seed tables as input to the WMLTP method to get the rehtted 
frequencies, i^nT- 

As we have demonstrated in Section 14.21 the fixed-point iteration that is used for the 
compilation of the custom-made seed tables which are required as input to the WMLTP 
method converges on a solution rather rapidly. Therefore, it is to be expected that the 
rehtted frequencies, generated in the second sub-step deviate only slightly from the 

observed frequencies, In fact, in most of our applications we have found that ~ 

As a result, the fourth step of our outlier correction scheme may be skipped provided the 
preceding hxed-point iteration has been carried out carefully. In this case, the remaining 
hfth step of onr outlier correction scheme is performed by using 


7.1.5. Step 5 

In the hnal step of our outlier correction scheme, we npdate all of the modes that have 
been hagged as being outliers in steps 1 through 3. In some cases we only npdate the 
frequencies of those outliers and in other cases we npdate only the associated uncertainties, 
and in still other cases we update both the frequencies and their uncertainties at the same 
time. In addition, as we shall demonstrate in Section 17.1.61 we also encountered a few cases 
for which neither their frequency nor their uncertainty have been altered by the procedure. 

Before making the updates, we first check whether the refitted frequencies actually 
constitute an improvement over the originally observed frequencies, by calculating the 
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distances (5oid and (5new, respectively, given by 


«id = 

c. = «r -)'. 


(64) 

(65) 


using the emulated frequencies and respectively, as reference values. If < 5, 


old’ 


we replace the original observed frequency, , with the rehtted frequency, for that 

mode, i.e., 

= <z", (66) 


obs.new 


and we set for the associated uncertainty 


Aqy™" = max (A<,-, A^f. Ad'S?, A^J'^ |<5‘ - <,»1 /3. Kf “ "”1/3) ■ (67) 


Here, the hrst four terms on the right-hand side represent all of the various uncertainties 
available for the mode {n,l) at this point, while the last two terms make sure that the 
rehtted frequency is off by no more than 3 a from the curve that was ht to the unsealed 
frequency differences, Az/„ ;, as well as from the curve that was ht to the scaled frequency 
diherences, in the hrst and second step of the outlier correction scheme, respectively. 

If A <5oid; do not adjust the original observed frequency, z/®*)®, but only replace the 
associated uncertainty with 
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Here, as compared to equation (j67)) the term Az/“®" has been dropped, and in the last two 
terms z/jj®'^ has been replaced with z/®*)®. 


7.1.6. Summary of outlying cases 

By means of our outlier correction procedure we detected a total of 2,483 outliers in the 
entire set of 6,366 modes that we originally inverted (i.e., a total of 39 % of the htted modes 
were hagged as being outlying cases). For 1,156 of them we found improved frequencies, 
z/“®'^, by rehtting the corresponding modes, and by then replacing the original observed 
frequencies, nfff, with z/“()" whenever while for the remaining 1,327 outliers it 

turned out that the corresponding rehtted frequencies were not an improvement over the 
original observed frequencies, z/®*)®, because of and hence we left the original 

observed frequencies unchanged. 

In Figure [25] we illustrate the location of the outliers in the l-u plane that were detected 
by means of our outlier correction procedure. The 1,156 outliers for which we were able 
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to compute an improved refitted frequency are shown as the red diamonds, while the 1,327 
outliers for which our rehtted frequency was not an improvement over our original, observed 
frequency, are shown as the green diamonds. When we examined the spherical harmonic 
degree distribution of both sets of outliers, we found that both sets were primarily located 
at low and moderate degrees rather than at high degrees. In fact, we found that 75 % of the 
2,483 outliers were located at or below I = 350. 

Of the 1,156 cases for which we were able to compute an improved rehtted frequency, 
we found far more cases (e.g., 1,031) for which the rehtted frequency was lower than the 
original frequency. That left only 125 cases for which the rehtted frequency was higher than 
was the original frequency. Furthermore, for only 143 of the 1,156 cases were we able to hnd 
a normalized frequency adjustment that was greater than 3 a in absolute magnitude. 

By means of either equation fl67)) or fl6S]) we tried to hnd improved estimates for the 
frequency uncertainties of the 2,483 outliers. Such adjustment of the frequency uncertainties 
is in a sense equivalent to the allowance of systematic errors in the observed frequencies, 
and respectively. A summary of the statistics that resulted when we carefully 

evaluated the frequency uncertainties of all 2,483 outliers using equations and fl65D 
is given in Table [71 In the hrst column of Table [7] the eight terms are listed that could 
trigger the modihcation of the frequency uncertainty of an outlying case according to either 
equation fl67|) or fl68|) . In the second column of Table [7] we list the number of cases for which 
one of the terms included in the right-hand side of equation (167)) triggered the adjustment 
of the frequency uncertainty. Interestingly, we found that for 286 cases, 

as shown in the hrst row of column 2 in Table [71 and therefore we did not alter the original 
frequency uncertainty, of those cases. For the remaining 870 cases we increased the 

original frequency uncertainty because one of the terms listed below in the hrst column 

of Table [Tlturned out to be larger in magnitude than In the third column of Table [7) we 

list the number of cases that we evaluated using equation ([65]) . For 1,129 of these 1,327 cases, 
we ended up increasing the original frequency uncertainties; however, for the remaining 198 
cases the original frequency uncertainty was large enough that it did not need to be replaced 
with one of the other terms in equation ([68]) . In summary, considering all of the 2,483 outliers 
we found an average ratio of the hnal and original uncertainties of 1.974 ± 1.237. When we 
included only the 1,999 cases for which the frequency uncertainty actually increased in the 
computation of the average ratio of the new to the original uncertainties, we found that the 
average ratio was increased slightly to 2.210 ± 1.271. 

The 198 outliers listed at the top of the third column in Table [7] are those cases for which 
neither the frequencies nor the frequency uncertainties needed to be changed. Therefore, the 
question arises as to why those 198 cases were flagged as being outliers in the hrst case. This 
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can happen because we are invoking two criteria for the identihcation of outliers which do 
not depend in any way upon the frequency uncertainty, These criteria are given in 

equation fl5B]) and in the third branch of equation respectively. As a result, if a mode is 
initially flagged as being an outlier by one of those two criteria, but later its rehtted frequency 
is not found to be an improvement over the original frequency of that mode, , and also 
if the original frequency uncertainty of that mode, is large enough in magnitude that 

equation (l68l) sets this mode ends up with both its original frequency 

and its original uncertainty intact. We note, however, that such a mode is not expected to 
hamper the convergence of the inversion code because of its exceptionally large frequency 
uncertainty, that has been returned by the WMLTP code for this mode. 

In spite of the above discussion about the number of outliers that our procedure identi- 
hed, we wish to stress that, in the majority of the 6,366 original hts for which v < 4500 /iHz, 
we did not hnd it necessary to alter our original htted frequencies and/or uncertainties 
thereof at all. We illustrate this important point further in Table [8] in which we have sum¬ 
marized the various outcomes of our hve-step procedure. In the second row of Table [8] we 
note that we did not have to alter the frequencies of 5,210, or 81.8 %, of our 6,366 original hts. 
This overwhelming majority of our cases was comprised of the 3,883 cases that were never 
identihed as being outliers by our procedure, the 1,129 cases for which only the frequency 
uncertainties needed to be altered, and the 198 cases for which neither the frequencies nor 
their uncertainties needed to be adjusted. Also, in the third row of Table Owe note that a 
total of 4,367, or 68.6 %, of our original cases did not need to have their uncertainties altered. 
This number was comprised of the 3,883 cases that were never identihed as being outliers, 
the 286 cases for which only the frequencies needed to be altered, and the 198 cases for which 
neither the frequencies nor the uncertainties needed to be altered. In the remaining rows 
of Table [8] we summarize the 2,483 outlying cases and their outcomes that resulted from 
the various steps of our procedure in diminishing order as percentages of our total of 6,366 
original cases. 

While Table [8] shows that our outlier identihcation and correction procedure ended up 
requiring an increase in a total of 1,999 of the original frequency uncertainties, we note that 
these increases were found to be needed because the original uncertainties that the WMLTP 
code computed were too small. In this regard, it is important to note that the formal 
uncertainties delivered by the WMLTP code only represent the statistical errors which are 
not corrected for any systematic ehects. Therefore, our outlier correction scheme is useful 
not only for detecting systematic errors in the sets of htted frequencies but also for detecting 
frequency uncertainties that are unduly small. We hnally note that it is both the correction 
of the htted frequencies for systematic errors as well as the adequate adjustment of the 
frequency uncertainties by means of our outlier correction scheme that hnally allowed the 
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inversion shown in Fignre[26]to converge properly on a solntion. 


7.2. Inversion results 


Accurate measurements of the high-degree solar oscillation modes are particularly im¬ 
portant for helioseismic diagnostics of the near-surface layers of the Sun. These layers are 
believed to play a critical role in the formation of the magnetic network, sunspots and active 
regions, and thus are a key to understanding the mechanisms of solar activity and variabil¬ 
ity. To illustrate the potential of the new method of measurements of the high-degree mode 
frequencies we performed a standard inversion procedure for the solar sound-speed prohle, 
which is based on the mode sensitivity kernels derived from the variational principle, and 
the OLA inversion method. The sensitivity kernels are calculated from the adiabatic eigen¬ 
functions, and the surface non-adiabatic effects are removed by assuming that these can be 
described as a function of frequency, weighted with the mode inertia (cf. equation (01])). 
The inversion method provides locally averaged estimates of the sound-speed variations at 
various target positions along the solar radius. The central locations and characteristic width 
(“spread”) of the averaging kernels of these estimates provide a quantitative measure of the 
resolving power of the obs erved frequency se t. The mathematical details of the inversion 
procedure are described bv iKosovichevI (119991) . 


The inversion results illustrated in Figure [26] show t he relative deviations of the square d 
sound speed in the Sun from the standard solar Model S fjChristensen-Dalsgaard et al.lll9961) . 
The results reproduce a well-known peak at the base of the convection zone, which is as¬ 
sociated with the tacho cline - a narrow area of the strong rotat ional shear at the bottom 
of the convection zone (jKosovichevIIl996al: lElliott fc Goughlll999l) . In the lower part of the 
convection zone, the sound-speed deviation is almost zero due to the adiabatic stratihcation 
of the radial structure. The sound-speed prohle gradually decreases relative to the model in 
the bulk of the convection zone, which may indicate a deviation from the adiabatic stratih¬ 
cation of the model (prescribed by the mixing-length theory of convection), and then shows 
a sharp decrease in the upper 5% of the solar radius (~ 35 Mm). Such decrease in the 
sound speed was previously detected by inversions, however, its structure was not resolved. 
The new results show that the structure of the near-surface layer is now well-resolved. The 
substantial deviation from the model coincide s with the subsurface rotational shear layer 
(called “leptocline” bv iGodier fc Rozelotll200ll) . and indicates that the structure of the solar 
convection zone and the heat transport properties may be signihcantly different from the 
predictions of the mixing-length theory. Of course, this result is only a starting point of 
systematic investigations of the global structure and dynamics of this layer, and also of the 


















variations with the solar cycle. 


8. Concluding remarks 

Global helioseismology has proven to be an extremely powerful tool for the investiga¬ 
tion of the internal structure and dynamical motions of the Sun and, thus, should not be 
abandoned in favor of local helioseismology alone. Numerical inversions of low- and medium- 
degree solar oscillation frequencies have conhrmed that the solar structure is in remarkably 
good agreement with predictions of the standard solar model. However, global helioseismic 
studies have not yet provided reliable information about the very interesting near-surface 
region of the Sun because high-degree modes were not measured accurately enough due to 
serious problems in the understanding of the structure of a ridge of oscillatory power and 
various distortion effects of the measuring device. 

In order to provide a more-accurate characterization of the high-degree ridges of power, 
we have presented in this paper a sophisticated mathematical method for the htting of m- 
averaged power spectra which employs a theoretical prohle that contains multiple peaks 
which represent each mode of interest and also its surrounding temporal and spatial side- 
lobes. This method includes both spatial leakage matrices which weight each of the spatial 
sidelobes differently according to the details of the instrument which provided the data from 
which the power spectra were created, and it also includes the temporal window function 
of each observational time series to allow for the data gaps that are usually present in the 
observations. 

In Section |2] above, we outlined several different approaches to the calculation of m- 
averaged spectra, and we also described the need for the correction of the frequency splitting 
coefficients which are used in the generation of such spectra for the effects of modal distortions 
that are introduced by solar latitudinal differential rotation. In Section [3] we presented 
mathematical details of the computation of the spatial leakage matrices which are required 
for the accurate htting of both isolated modal peaks and ridges of power, and we presented 
the mathematical details of the method by which such leakage matrices are corrected for the 
effects of differential rotation. In Section 0] we presented the mathematical details of this 
new method, which we have named the WMLTP method, and which we have also referred 
to in the text as our Method 2. 

In Section [5] we presented the results of a series of studies of the sensitivity of the £t- 
ted frequencies that this method produces to various details of the method of generating 
the m-averaged power spectra and of the approximation that is made to the shape of the 
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various leakage matrix peaks. In Section 15.1.11 we showed that the weighting of the indi¬ 
vidual tesseral and sectoral power spectra at each degree prior to the computation of the 
m-averaged spectrum had only a very minor influence on the resulting frequencies. On the 
other hand, in Section 15.1.21 we showed that the use of frequency-splitting coefficients that 
were computed in a non-n-averaged manner and which were also corrected for the distorting 
effects of differential latitudinal rotation dehnitely did result in substantial changes to the 
htted frequencies. In Section [5.1.3l we showed that the majority of these systematic frequency 
changes were due to the use of such non-n-averaged frequency-splitting coefficients in place 
of coefficients that were computed in an n-averaged manner. That is, we demonstrated the 
importance of using a narrow frequency range that was centered on each individual ridge 
when cross-correlating the non-zonal power spectra with the zonal spectrum at each degree. 
In Section ESI we showed that an increase of 18 % in the width of the Gaussian function that 
we employ in Method 2 as an approximation of the peaks in the effective leakage matrices 
did not result in any systematic frequency changes. 

In Section [6] we demonstrated the substantial improvements that the WMLTP method 
makes in the frequencies, linewidths, amplitudes, and in the uncertainties of all of these 
quantities in comparison with the same quantities computed using our earlier Method 1, 
which only employed a single Lorentizan prohle in place of the more complicated multiple- 
peak prohle of Method 2. In particular, we have pointed out that the use of a multiple-peak 
prohle in the WMLTP method reduces the formal frequency uncertainties by a factor of 15.2± 
0.3 on average, and it also reduced the formal linewidth uncertainties by a factor of 29.8±0.9 
on average. We have also demonstrated that the frequencies, linewidths, and amplitudes that 
we computed using the multiple-peak prohle do not exhibit any sharp discontinuities along 
the various ridges as did the same quantities that were obtained using the single-peak prohle. 
These discontinuities occurred where the individual peaks became blended into broad ridges 
of power at the higher degrees. We have also compared the relative smoothness of both the 
frequencies and linewidths which were computed with the WMLTP method from various 
sets of power spectra that were obtained with the MDI instrument. We have also displayed 
the sets of 12,359 frequencies, linewidths, amplitudes, asymmetries, and power levels that 
we generated by applying the WMLTP method to the 66-day long 2010 MDI Dynamics Run 
power spectra. 

In Section ITT] we presented the details of a new procedure that we developed to improve 
both the numerical stability and reliability of the structural inversion that we presented in 
Section 17.21 and of similar structural inversions that we expect to compute with similar 
Method 2 frequency tables in the future. This procedure both identihes and corrects outliers 
in both the frequencies and their associated uncertainties. We applied it to the subset of 
6,366 Method 2 frequencies that we had originally inverted and we obtained a much smoother 
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inverted sound speed deviation profile in the process. 


In conclusion, we have shown that the development of the WMLTP method of power 
spectral htting provides a powerful new tool for the provision of accurate and reliable estima¬ 
tion of low-, intermediate, and high-degree mode parameters. We demonstrated this claim 
impressively in Section 17.21 with the presentation of a new structural inversion correspond¬ 
ing to the beginning of Solar Cycle 24 that we computed using our cleaned table of MDI 
2010 Dynamics Run frequencies. This new inversion resolves for the hrst time the seismic 
properties of the upper convective boundary layer, and shows a substantial and surprisingly 
sharp deviation from the ad i abatic sound-speed prohle of the Standard Solar Model S of 
Christensen-Dalsgaard et al.l fll996l) . which is likely due to the near-surface turbulence ef¬ 
fects. Moreover, this is also evidence that the predictions of the mixing-length theory may 
be signihcantly different from both the actual structure of the convection zone and the heat 
transport properties in the subsurface layers of the Sun. By comparing upcoming structural 
inversions similar to that presented here in Section 17.21 with realistic 3D MHD numerical 
simulations of the upper solar convection zone and the subsurface shear layer we hope to 
make signihcant progress in the understanding of the structure and dynamics of the outer 
layers of the Sun. In this regard, we note that the detailed measurements of the linewidths 
of the /- and p -modes could also provide useful constraints on the properties of near-surface 
convection fcf. IChristensen-Dalsgaard et al.lll989l: iRabello-Soares et al.lll999l) . 


Because for degrees I ^ 1 m-averaging of the originally generated zonal, sectoral, and 
tesseral power spectra signihcantly improves the signal-to-noise ratio, the WMLTP method 
is particularly well suited for the htting of m-averaged power spectra derived from time series 


shifts with the phase of t he solar activity cycle (see, e.g., Ronan et ah Il994 : Jeheries 


Rhodes et al 


Jain et al. 


2002 


2011 


2003 


iRose et al. 


20031: iRabello- Soares et al 


1994; 

Jeheries 

1998 

Trinathv et al. 

2010 


Rhodes et al.ll201l|) . On the other hand, the generation of m-averaged 


power spectra requires knowledge of suitable frequency splitting coefficients (cf. Sect. 12.2p . 
Therefore, as we mentioned in Section 12.31 we have also begun to develop our MPTS htting 
method which operates directly upon all 2/ -|- 1 un-averaged power spectra for each degree. 
The MPTS method allows us to determine all 2/ -I- 1 sets of modal parameters of a given 
(n, 1) multiplet. Once we have determined all 2Z -|- 1 of the frequencies of that multiplet, 
we can then compute both the average frequency and the frequency-splitting coefficients 
for that multiplet. In essence, the MPTS method will allow us to generate sets of non-n- 
averaged frequency splitting coefficients without recourse to the cross-correlation method 
that we described in Section [221 Details of the MPTS method will be given in the upcoming 
second paper of our series of three papers. 


























































Unfortunately, as we mentioned earlier the MPTS method cannot be applied to power 
spectra that are computed from time series that are as short as three days in duration, nor 
can it £t peaks having frequencies greater than about 5000 /iHz because of the low signal- 
to-noise ratios of the underlying power spectra. As a result of these two limitations, both 
the WMLTP and MPTS htting methods are complementary to each other. Hence, we hope 
that, once we have been able to bring the current version of the MPTS code up to the level 
of the rev6 version of the WMLTP code, the parallel application of both htting methods will 
lead to new information concerning temporal variations in the Sun’s sub-surface structure 
and dynamics. 


In this work we utilized data from the Solar Oscillations Investigation / Michelson Doppler 
Imager (SOI/MDI) on board the Solar and Heliospheric Observatory (SOHO), and we have 
made use of NASAs Astrophysics Data System. SOHO is a project of international cooper¬ 
ation between ESA and NASA. The SOI/MDI project is supported by NASA grant NAG5- 
10483 to Stanford University. The portion of the research that was conducted at the Uni¬ 
versity of Southern California was supported in part by NASA Grants NNX08A24G, NAG5- 
13510, NAG5-11582, NAG5-11001, NAG5-8545, NAG5-8021, NAG5-6104, and NAGW-13, 
by Stanford University sub-awards 14405890-126967, 1503169-33789-A, and 29056-C, by 
Stanford University Sub-Contract Number 6914, and by USC’s Office of Undergraduate Pro¬ 
grams. Part of this work is the result of research performed at the Jet Propulsion Laboratory 
of the California Institute of Technology under a contract with the National Aeronautics and 
Space Administration. We thank the anonymous referee for his valuable contributions to im¬ 
prove the presentation of this work. J.R. is grateful to R. Bulirsch, P. Rentrop, and B. Vexler 
of the Technische Universitat Miinchen for their generous support and hospitality, and to 
K. Schittkowski of the University of Bayreuth for providing the source code of his NLPQL 
optimization technique. 


A. The single-peak, averaged-spectrum method 


The Single-Peak, Averaged-Spectrum, or SPAS, Method, which is also referred to as 
Method 1, was our initial, hrst-generation htting method. It was developed in the late-1980s 
through the mid-1990s in order to ht the peaks in the MDI power spect ra that were generated 
as part of that experiment’s Structure Program fjScherrer et al.lll995l) . This work led to the 
publication of a set of tables of /- and p-mode frequencies and their associated errors which 
has become a major reference in the helioseismic literature for the MDI Medium-/ Program 
flRhodes et al.lll997l) . We also have applied Method 1 to the a nalysis of power spe ctra which 
were derived from MDI’s Full-Disk (FD) Dynamics Program fjScherrer et al.lll995l) . This led 
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to the publication of 
the MDI FD spectra fiRhodes et alJll998alJbl) . 


he hrst frequencies to be computed for the high-degree modes from 


In Method 1 a single, symmetric Lorentzian prohle plus a linear background term is 
used as the model prohle to represent an oscillation peak, i.e.. 


M{u,p) = 


A 


+ a + bi/, X = 


2 (z/ - z/q) 


1 + 

rpi 

p= {A,uo,w,a,b) , 


w 


(Al) 

(A 2 ) 


with z/ being frequency. The hve ht parameters are collected in the vector p, and include, 
respectively, the mode amplitude A, the mode frequency z/q, the mode linewidth w, and 
the background noise parameters a and b. The vector p is determined by htting the model 
prohle flAlIl in the least-squares sense to the m-averaged spectrum in a htting range centered 
about the peak of interest that is selected by using the heuristic approach as described below. 
The resulting unconstrained least-sq uares problem is s olved by using the FORTRAN routine 
Imderl from the MINPACK project flMore et al.lll980l) . This routine is designed to minimize 
the sum of squares of Ai nonlinear fu nctions in A? va r iables fA? < AO by a modihcation of 
the Levenberg-Marquardt algorithm fjLevenberdll944J: lMarquardtlll963l) . 


The variance of the ith ht parameter p, G p is calculated from 

2F 


varpi = 


Al — A2 


H, 


-1 


(A3) 


where Ai is the number of data points (frequency bins), A 2 is the number of ht parameters. 


is the ith diagonal elemen t of the inverse of the Hessian matrix if, and S is the sum 


of squares of the ht residuals flNAG Fortran Library. Mark 231120111) . Both H and S are 
calculated at the solution of the least-squares problem. In the neighborhood of the solution 
the Hessian matrix can be adequately approximated by 


H = 2 rj, 


(A4) 


with J being the Jacobian matrix, thereb y avoiding the n eed to compute or approximate sec¬ 
ond derivatives of the objective function flGill et al.lll98ll) . The 100 (1 — 7 ) percent conhdence 
interval on the fth ht parameter pi is given by 


Pi ± Vvarpi fi_.^/2,Ai-A2, 


(A5) 


where fi_.y/ 2 ,Ai-A 2 is the 100 (1— 7 / 2 ) percentage point of Student’s (-dist ribution with Ai —A 2 
degrees of freedom flWolberdll967l: INAG Fortran Library. Mark 23ll201ll) . For 1 a conhdence 
intervals 7 = 0.31731, and for Ai — A 2 ^ 1 we have fi_.y/ 2 ,Ai-A 2 ~ 1- 
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We have implemented into Method 1 an option which allows to use either a narrow 
fitting range or else a wide fitting range. We switch between the narrow and the wide fitting 
range based upon the ratio, r, between the linewidth and the frequency change with respect 
to degree, namely 


w 

Au/Al 


(A6) 


We estimate this ratio using linewidths and frequencies from seed tables that have been 
derived from previously computed modal parameters. We use the narrow fitting range when 
< ■'"crit and the wide htting range for r > rent, namely when the modes are well separated, 
or not, from the spatial leaks of nearby degrees. Because there is no sharp transition between 
modes and ridges of power, the choice of rent is somewhat arbitrary. By visual inspection of 
the hts obtained with Method 1 we found the values 


0.65 for n = 0, 
0.80 for n > 0, 


(A7) 


to be useful. For the n = 0 ridge we had to select a separate value of rent because for this 
ridge individual modal peaks blend into ridges of power for smaller ratios, r, than for the 
other ridges. The values of rerit given in equation flA7p correspond to I ~ 278 for u = 0, 
I 208 for n = 1, and / 15 for n = 29. 

In the left panels of Figure [14] we show typical hts obtained with Method 1 when applied 
to the m-averaged spectral set iS2010_66a. The hts obtained by using a narrow htting range 
are shown as the green lines, while the blue lines are for hts that employed the wide htting 
range. We note that the use of both the narrow and wide htting ranges for [n, 1) = (2, 70) 
and (2, 200) in the top-left and middle-left panel of this Figure required the cancellation of 
the criterion (]A7p . which would have enforced the use of the narrow htting range for (2, 70) 
and the use of the wide htting range for (2, 200). From the colored tick marks in these two 
panels it becomes evident that the frequencies determined with the wide htting range diher 
signihcantly from those obtained with the narrow htting range. This is also demonstrated 
in Table [6l We stress that while this disagreement is obvious in the cases of those modes 
for which the peaks are well resolved, the very same problems are occurring in the cases of 
broad ridges of power even though the observed peaks do not show evidence of the individual 
spatial leaks as is shown here in the bottom-left panel of Figure [TH 


This so-called “frequency pulling” of the measured ridge-ht frequencie s away from the 
true solar frequencies was hrst addressed by iLibbrecht fc Kaufman! fll988l) , and is mainly 
caused by the asymmetry of the power distribut ion of the spatial leaks with respect to 
the target mode frequency flKorzennik et al.ll2004l) . We first tried to resolve this problem 
by developing a frequency-correction scheme in which we htted those modes for which the 
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target peak and the spatial leaks are well separated with both the narrow htting range and 
the wide htting range in Method 1 as is shown here in the top-left and the middle-left panel 
of Figure [T3J In order to be able to do so we had to override criterion flA7p . of course. 
By means of a multiple linear regression we then tried to correct the frequencies obtained 
with the wide htting range for the ehects of mode blending using the diherences between 
the narrow- and wide-ht frequencies. However, we eventually had to abandon this approach 
because the range of degrees for which the multiple linear regression model can be directly 
computed, namely, those degrees for which the modal peaks are not blended into ridges of 
power, is not large enough to allow the regre ssion model to be safely extrapolated up to 
degrees of 1000 and higher where it is needed flRhodes et al.ll200ll) . 


B. Estimates of smoothness of fitted parameters 

On physical grounds it can be assumed that any modal parameter, e.g., frequency or 
linewidth, is a smooth function of spherical harmonic degree / along a ridge of radial order 
n. We suggest to measure the smoothness of any modal parameter along a ridge by the 
so-called normalized point-to-point scatter, S{v), dehned by 

-I ^max 1 

“"f") = TT, — zrr-\ E - "‘t. (Bi) 

^ \^m&x ^m'm) 

^—^min 

where v is the variable/parameter whose smoothness is estimated, and /min and /max are, 
respectively, the minimum and maximum degree / of the portion of interest of the ridge. 
Because S{v) is rather sensitive to missing cases Vi we recommend that any gaps that might 
exist in the range from /min to /max are hlled by interpolation. We note that the value of 
Z'(t;) is biased towards any slope of the variable v that might exist in the range from /min to 
/max because such slope causes the differences Uj+i — Vi in equation flBlD to be systematically 
different from zero. In general, however, this bias is concealed by the scatter of the differences 
Uj+i — Vi due to noise. 
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Table 1. Epochs of observing runs used in this work. 


run 

starting and end date 

duration 

(days) 

Duration 
(60-sec samples) 

gap-hlled 
duty cycle 

7^1996_61 

May 24 - Jul 23, 1996 

60.75 

87,480 

97.30 % 

7^2001_90 

Feb 24 - May 28, 2001 

90 

129,600 

97.02% 

7^2010_66 

May 7 - Jul 11, 2010 

66 

95,040 

93.87% 

7^2010_03 

May 7 - May 9, 2010 

3 

4,320 

85.52% 


Note. — Column 1 lists the naming convention we assigned to each observing 
run. Columns 2 through 5 contain the starting and ending dates, the durations 
in days, the durations in 60-second samples, and the gap-hlled duty cycles of the 
four runs, respectiv ely. Each observing run was obtained as part of the MDI 
Full-Disk Program flScherrer et al.l 119951) . 
















Table 2. Sets of m-averaged spectra used in this work. 


observing 

run 

m-averaged 
spectral set 

m-averaging 

frequency splitting 
coefficients 

7^1996_61 

51996_61 

unweighted 

corrected, n-averaged 

7^2001_90 

52001_90 

unweighted 

corrected, n-averaged 

7^2010_66 

52010_66a 

unweighted 

corrected, n-averaged 

7^2010_66 

52010_66b 

weighted 

corrected, n-averaged 

7^2010_66 

52010_66c 

weighted 

raw, non-n-averaged 

7^2010_66 

52010_66d 

weighted 

corrected, non-n-averaged 

7^2010_03 

52010_03 

unweighted 

corrected, n-averaged 


Note. — Column 1 lists the observing run using the naming conven¬ 
tion introduced in Table [U Column 2 lists the naming convention we 
assigned to each m-averaged spectral set that we computed from the 
un-averaged power spectra obtained from the observing run listed in 
the same row of column 1. Column 3 indicates whether the averaged 
spectra were computed using the weights described in Section 12.21 or 
were computed in an unweighted manner. Column 4 describes the 
type of frequency splitting coefficients that were employed in the av¬ 
eraging step. 
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Table 3. Performance of the fixed-point iteration for the construction of accurate seed 
tables using the example of the frequency, v. 



Note. — Columns 2 through 5 contain the number of cases considered, 
Tic, the number of rejected outliers, Uout, the average, avg, and the standard 
deviation, std, of the cases listed in column 1. The values of the frequency 
derivative, /c = 1, 2, 3, were input as seeds into the WMLTP 

code to generate the fitted frequencies, ^ ^ ~ 15 2,3. The differences 

of the frequency derivatives, (Au /that are listed in rows 1 and 
2, respectively, are unsealed, while the differences of the fitted frequencies, 
that are listed in rows 3 and 4, respectively, are scaled by the maxi¬ 
mum of the uncertainties of the frequencies that are subtracted. 
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Table 4. Fitted parameters invoked in the individnal passes used for determining the 

htting vector p. 


pass 

htted parameters 
main peak & background 

htted parameters 
n-leaks 

#1 

^n,l 



a, b, 

c 




#2 


l^n,l 

'^71,1 






#3 

An,l 

Vn,l 

'^n,l 

a, b, 

c 




#4 

An,l 



a, b, 

c 

A. 



#5 


l^n,l 


a, b, 

c 




#6 



'^n,l 

a, b, 

c 



Wi 

#r 

^n,l 

^n,l 


a, b, 

c 

A^ 

Vi 

Wi 

#8 

^n,l 

Vn,l 


a, b, 

c 

Ai 

Vi 

Wi 

#9 

An,l 

l^n,l 

'^n,l 

a, b, 

c 

A, 

Vi 

Wi 

#10 

An,l 

l^n,l 

'^n,l 

a, b, 

c 

A, 

Vi 

Wi 

#11 




Bn,l 




B, 

#12 

An,l 

Vn,l 

'^n,l 

a, b, 

c 

A. 

Vi 

Wi 

#13 

An,l 

Vn,l 

'^n,l 


c 

A. 

Vi 

Wi Bi 


Note. — The value of a parameter that is not invoked 
in a specihc pass is set to either its initial guess or else to 
its value as determined in a previous pass. For the initial 
guess of the line asymmetry parameter we use = 0 and 
Bi = 0, i = 1,..., N. The initial guesses of the remaining 
£t parameters are taken from a seed table. If the symmet¬ 
rical prohle is used Bn,i and Bi, i = 1,..., N are kept hxed 
to zero in passes #11 and #13. To increase the numeri¬ 
cal stability of the WMLTP code rather tight constraints 
are applied to the htted parameters in passes # 1 through 
# 7. The rationale for passes # 8 through passes # 10 is to 
gradually relax these constraints as much as possible. 
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Table 5. Statistical analysis of the magnitude of raw and scaled frequency differences 
obtained from four comparisons Cl through C4 in the study of the sensitivity of Method 2 
to the details of the m-averaging procedure and the effective leakage matrix, respectively. 




u < 7000 /iHz 

V < 4500 /iHz 



avg ± std 

ste 

"^tot 

"^out 

avg ± std 

ste 

^tot 

^out 


|Az/| 

0.265 ±0.514 

0.005 

12,353 

55 

0.044 ±0.137 

0.002 

6,421 

4 

|Az//(TAi2| 

1.062 ± 1.075 

0.010 

12,353 

3 

1.076 ± 1.177 

0.015 

6,421 

3 


|Az/| 

0.403 ± 0.484 

0.004 

12,158 

24 

0.278 ±0.301 

0.004 

6,235 

0 

|Az//aAi.| 

4.858 ±6.398 

0.058 

12,158 

79 

7.786 ± 7.689 

0.098 

6,235 

79 


|Az/| 

0.328 ± 0.403 

0.004 

12,159 

38 

0.302 ±0.331 

0.004 

6,236 

0 

|Az//(TAi.| 

5.786 ± 8.403 

0.076 

12,159 

5 

10.225 ±9.759 

0.124 

6,236 

5 


|Az/| 

0.149 ±0.352 

0.003 

12,368 

38 

0.060 ±0.087 

0.001 

6,411 

2 

|Az//(TAi.| 

1.000 ± 1.247 

0.011 

12,368 

2 

1.370 ± 1.452 

0.018 

6,411 

2 


Note. — Cl denotes the comparison of the frequencies from table J^2010_66a with the corre¬ 
sponding frequencies from table J^2010_66b; C2 denotes the comparison of the frequencies from 
table J^2010_66c with the corresponding frequencies from table J^2010_66d; C3 denotes the com¬ 
parison of the frequencies from table J^2010_66b with the corresponding frequencies from table 
J^2010_66d; and C4 denotes the comparison of the frequencies from table J^2010_66b with the cor¬ 
responding frequencies from table J^2010_66b.lkm. For each of the four comparisons Cl through 
C4 we show the average, avg, the standard deviation, std, and the standard error of the average, 
ste, of the magnitude of the raw frequency differences, |Az/|, measured in pHz, in the hrst row, and 
of the magnitude of the normalized frequency differences, |Ai//(TA!y|, where ctai/ denotes the formal 
error of the frequency difference Ai^, in the second row. In our analysis we weeded out from the 
given ritot frequency differences Uout gross outliers by using the rejection criterion |Az/| > 5/iHz 
for the raw differences and |Az//crAj/| > 35 for the normalized differences. As shown in the last 
four columns, we repeated the entire analysis for the subset of modes having v < 4500 /iHz. 
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Table 6. Normalized comparison of sample Method 1 and Method 2 £t results based upon 
the m-averaged spectral set iS2010_66a for the modes (n, 1) = (2, 70), (2, 200), and (2, 600). 


(n,/) 

wfr ~ i^l,nfr 

^2 - 1^1,nfr 

0'u,2 

1^2 — i^l,wfr 

O'ua 

'W^l,wfr — "U^l.nfr 

W2 - Wi,nfr 

^w,2 

W2 — Wi^wfr 

^w,2 

(2,70) 

862.7 

0.2 

-258.6 

2225.0 

24.9 

-1121.3 

(2,200) 

(2,600) 

10.1 

-10.6 

-210.8 

-16.4 

10.6 

5.3 

-669.3 

-261.8 


Note. — Column 2 contains the normalized frequency differences between the narrow (nfr) 
and wide (wfr) htting range versions of Method 1 for which the denominator was narrow 
htting range uncertainty. Column 3 contains the normalized frequency differences between 
Method 2 and the narrow htting range version of Method 1 for which the denominator was 
the Method 2 uncertainty. Column 4 contains the normalized frequency differences between 
Method 2 and the wide htting range version of Method 1 for which the Method 2 uncertainty 
was the denominator. Columns 5 through 7 contain the similar quantities for the linewidths. 
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Table 7. Classification of the 2,483 outliers detected by means of the outlier correction 
scheme, as described in Section 17.11 in terms of the various terms that could trigger the 
adjustment of the frequency uncertainty of the outlying cases. 



no. of cases 

term 

eq. dnZD 

eq. (ESD 

Ai'S" 

286 

198 


103 



13 

5 


224 

209 

Kf - <ri/3 

329 


1"™' - <i'l/3 
Kf - <7\/3 

201 

482 

Kf - <,'l/3 


433 

total 

1,156 

1,327 


Note. — In the first column we 
list the terms as given in the right- 
hand side of equations (j67|) and ([68]), 
respectively, that trigger the evalu¬ 
ation of the frequency uncertainties. 
The meaning of the various terms is 
given in Section 17.11 In column 2 
we list the number of cases that were 
triggered by the corresponding term 
listed in column 1 in case of equa¬ 
tion fl67|) . We note that this equa¬ 
tion is for those cases for which we 
were able to compute an improved re¬ 
fitted frequency in the manner as de¬ 
scribed in Section 17.1.41 Column 3 
contains the same quantities as col¬ 
umn 2, but the quantities in this col- 
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umn were triggered by the different 
terms in the right-hand side of equa¬ 
tion fl68l) . which applies for those cases 
for which our refitted frequency was 
not an improvement over our original, 
observed frequency. In the last row 
we give the total number of cases that 
were triggered by the respective equa¬ 
tions. A horizontal bar in a column 
means that the corresponding term 
listed in the first column does not ap¬ 
ply for the respective equations. 
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Table 8. Statistical overview of the results of the application of our 5-step outlier 
identihcation and correction procedure to the table of hts that was employed in the 

structural inversion shown in Figure [221 


# 

% 

Description 

6,366 

100.0 

Original hts for which v < 4500 /rHz 

5,210 

81.8 

Cases for which v was unchanged 

4,367 

68.6 

Cases for which Az/ was unchanged 

3,883 

61.0 

Cases that were never identihed as being outliers 

2,483 

39.0 

Cases that were identihed as being outliers 

1,999 

31.4 

Outliers that resulted in an increase in Az/ 

1,327 

20.8 

Outliers for which v was unchanged 

1,156 

18.2 

Outliers for which v was altered 

1,129 

17.7 

Outliers for which v was unchanged, but Az/ was increased 

1,031 

16.2 

Outliers for which v was decreased 

870 

13.7 

Outliers for which both v and Az/ were altered 

286 

4.5 

Outliers for which v was altered but Az/ was unchanged 

198 

3.1 

Outliers for which neither v nor Az/ were altered 

125 

2.0 

Outliers for which v was increased 


Note. — Here, v denotes the frequency and the uncertainty thereof. 
All of the percentages shown in the second column were computed using the 
6,366 original cases in this table. 
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Fig. 1.— Location rt of the inner turning point (a), and depth of penetration R — rt (b), in 
units of the solar radius R, for p-modes in a standard solar model. 
as fun ctions of degree I, for three typical frequencies. Adapted from 

(I2nn3h . 
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Fig. 2.— (left) Degree dependence of the n-averaged odd-order splitting coefficients oi, 03 , 
and 05 , in the range 5 < / < 1000 computed from the set of un-averaged power spectra 
obtained from the 771996_61 observing run, with the Oi coefficient at the top. As can be 
clearly seen, these raw splitting coefficients illustrate large discontinuities at about I ~ 210 . 
(right) Same as shown in the left panels, but after the discontinuities have been removed 
using the adjustment procedure described in Section \2?2[ 
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Fig. 3.— Dependence of m-averaged leakage matrix for a displacement component ratio of 
Ct = 1.0 on the difference, Al = I' — I, between the degree /' of the spatial side-lobes and 
the degree / of the target mode. The left panel is for (n,/) = (0,400), while the right panel 
is for (2,1000). The diamonds represent the original leakage matrix, while the triangles are 
for the leakage matrix which is corrected for the effect of the solar differential rotation. The 
full black lines are Gaussian hts to the uncorrected leakage matrices, and the dashed black 
lines are Gaussian hts to the corrected leakage matrices. The width of the corrected matrices 
increases considerably with increasing degree 1. In each panel the dashed red line represents 
a Gaussian the full-width-at-half maximum (FWHM) of which is larger by a factor of 1.18 
than the FWHM of the Gaussian shown by the black dashed line. These such distorted 
Gaussians are used in Section [52] to study the effect of changes in the FWHM of the leakage 
matrix upon the htted mode parameters. 
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Fig. 4.— (left) Dependence of the parameter which according to equation fl20|) is 

related to the width of the leakage matrix, on the offset, Xm’^\ of the leakage matrix due to 
the horizontal component of the modal velocity eigenfunction for (n,/) = (11,10) (green), 
(0,400) (red), and (2,1000) (black). The diamonds are for hts to a numerical representation 
of these leakage matrices that were corrected for the solar latitudinal differential rotation. 
The solid lines are for the parameterization as given in equation fl22l) . We particularly note 
that for high degree values the parameter is nearly independent of the value of Xm'^\ 
(right) Same as left panel, but for the solar velocity eigenfunction component ratio, In 

this panel the solid lines are for the parameterization given in equation ([23]). 
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Fig. 5.— Flow chart of the individual steps performed in the hxed-point iteration for the 
construction of accurate seed tables using the example of the frequency, v. Similar flow 
charts apply to amplitude and linewidth, respectively. 
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Fig. 6.— Half of the width, Wb, of the htting box measured in terms of (Au/An) versus 
degree for a selected set of radial orders, n. Here, (Au/An) denotes the mean value of 
and Az/„^;/A?7,| right with Az/„^;/A77,|ieft and Az/„^//An| right being the variation of 
the target mode frequency with respect to the radial order, n, at the lower and the higher 
frequency side of the htting box, respectively. For most cases shown the htting box widens 
with increasing degree by only about 10 to 30%, while it widens by about 150% for the 
n = 0 ridge. However, for the n = 8 and n = 29 ridge the width decreases with degree by 
about 23 % and 15 %, respectively. For the n = 6 ridge the width is practically constant. 
For the n = 20 ridge the n-leaks (n ± 1, /) and (n ± 2, /) are located within the htting box 
for higher degrees. 
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Fig. 7.— (left) Line asymmetry parameter B versus frequency for the n = 6 ridge computed 
from the m-averaged spectral set iS2010_66a. The error bars shown for selected values of 
frequency are for 1 a errors. Because the excursions of B to both positive and negative values 
for u > 5000 /iHz are hard to explain on physical grounds, we switched over to a symmetric 
Lorentzian prohle for the high-frequency portion of each ridge, (right) Same as left panel, 
except that the symmetrical profile was used for frequencies u > 4815 /iHz for this ridge. In 
both panels the red dashed line is for B = 0. 
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Fig. 8.— Frequency dependence of the numerical derivative of the frequency with respect to 
degree, Au/Al (in /iHz), for the n = 4 ridge. In the left panel the asymmetrical prohle has 
been used for htting all of the cases for this ridge, while in the right panel the asymmetrical 
prohle has been replaced with the symmetrical prohle for frequencies u > 4700 /iHz. Both 
sets of hts were computed from the m-averaged spectral set iS2010_66a. The error bars that 
are shown for selected cases are 1 a errors. The dashed red line is for Au/Al = 0. 
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Fig. 9.— (left) Straight line fit to the spectral power density in the troughs of the spectrum 
of the (n, 1) = (24, 8) mode, using the example of the m-averaged spectral set iS2010_66a. 
The spectral power density is shown in black, while the htted line is shown in red. The slope 
of this htted line is used as an estimate of the slope, b, of the background portion of the 
theoretical model prohle given in equations (|2T)) through fl29l) . when this prohle is htted to 
the spectrum in the htting box indicated by the vertical dashed lines in magenta. The htted 
prohle is shown as the green line, and the vertical blue dashed line is for the resulting htted 
frequency, (right) Same as left panel, except that it is for the (24, 56) mode. 
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Fig. 10.— (upper-left) Frequency dependence of the frequency differences, Au, which re¬ 
sulted from frequency tables J^2010_66a and J^2010_66b, in the sense of “weighted” minus 
“unweighted”. In each case, the same set of corrected, n-averaged frequency splitting coef- 
hcients has been used for the collapsing of the spectra, (upper-right) Degree dependence of 
the set of frequency differences, Au, shown in the upper-left panel, (lower-left) Frequency 
dependence of the normalized frequency differences, Avjo/^y. The normalization was carried 
out by dividing the raw frequency differences, Ai/, as shown in the upper-left panel, by the 
formal error, of each difference, (lower-right) Dependence of the low-frequency subset of 
the normalized frequency differences upon the fractional inner turning-point radii, rt/RQ, of 
those modes. Here, Rq is the radius of the Sun. Only mode frequencies less than 4500 /iHz 
were included in this panel. In all four panels the red dashed line is for a frequency difference 
of zero. In the two panels in the lower row the green dashed lines show the ±3 a values. 
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Fig. 11.— Same as Figure [101 except that it shows the dependence of both the raw and 
normalized frequency differences, Ai/, which resulted from frequency tables J^2010_66c and 
J^2010_66d, in the sense of “corrected” minus “raw” with respect to frequency, degree, and 
fractional inner turning-point radius. 
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Fig. 12.— Same as Figure HOl except that it shows the dependence of both the raw and 
normalized frequency differences, Ai/, which resulted from frequency tables J^2010_66b and 
J^2010_66d, in the sense of “non-n-averaged” minus “n-averaged” with respect to frequency, 
degree, and fractional inner turning-point radius. 
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Fig. 13.— Same as Figure HOl except that it shows the dependence of both the raw and 
normalized frequency differences, Ai/, which resulted from frequency tables J^2010_66b and 
J^2010_66b.lkm, in the sense of “wider Gaussian” minus “original Gaussian” with respect to 
frequency, degree, and fractional inner turning-point radius. 
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Fig. 14.— Fits to segments of three different spectra from the m-averaged spectral set 
52010_66a that were centered around the frequency of the (u, 1) = (2, 70) mode (top panels), 
the (2, 200) mode (middle panels), and the (2, 600) mode (bottom panels). In each panel the 
black line is for the m-averaged spectrum. In the top-left and middle-left panel the green line 
is for the £t using the narrow fitting range in Method 1, while the blue line is for the fit using 
the wide htting range in Method 1 to simulate the effect of fitting broad ridges of power. 
The segment centered on the (2, 600) mode has been ht using only the wide fitting range in 
Method 1, and the fitted prohle is shown as the blue curve in the bottom-left panel. All of the 
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computed profiles shown in the left-hand panels were generated using the symmetric prohle. 
The red lines in the th ree right-hand panels w ^ere al l computed using Method 2 employing the 
asymmetric prohle of iNigam fc KosovichevI (1l998h . For all six panels the htted frequencies 
are indicated by the colored tick marks that are located along both the top and bottom axes 
of each plot. 
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Fig. 15.— (upper-left) Frequency dependence of the ratios of the Method 1 narrow htting 
range frequency uncertainties divided by the Method 2 frequency uncertainties, (upper-right) 
Frequency dependence of the ratios of the Method 1 wide htting range frequency uncertainties 
divided by the Method 2 frequency uncertainties, (lower-left) Frequency dependence of the 
ratios of the Method 1 and Method 2 linewidth uncertainties for the narrow htting range 
cases, (lower-right) Frequency dependence of the ratios of the Method 1 and Method 2 
linewidth uncertainties for the wide htting range cases. In all four panels the dashed green 
lines show the average ratios, while the dashed red lines show the error ratios of unity. Note 
that the vertical scales are diherent in all four panels. In all four panels the m-averaged 
spectral set 52010_66a has been htted. 
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Fig. 16.— (upper-left) Frequency dependence of two sets of frequency differences, Au = 
j/obs _ j^mod^ between our fitted frequencies, and the corresponding Model S frequencies, 
j^mod, n = 0 ridge. The Au values as computed by Method 1 using the narrow 

htting range are shown as the black line, while the Au values computed using Method 2 
are shown as the red line. The outlier in both the Method 1 and Method 2 frequencies at 
about 1500 /iHz is for I = 219, and is caused by a glitch in the underlying spectrum. In 
both upper panels the dashed blue line represents a difference of zero, (upper-right) Same 
as upper-left panel, except that the wide htting range was used for computing the Method 1 
frequencies, (lower-left) Frequency dependence of the linewidths computed using Method 1 
with the narrow htting range (the black line) and using Method 2 (the red line). The dashed 
blue line represents linewidths of zero, (lower-right) Frequency dependence of the ratio of the 
linewidths as computed using Method 2 divided by the linewidths computed using Method 1 
with the wide htting range. The dashed line represents a linewidth ratio of unity. In all four 
panels, the hts were computed using the m-averaged spectral set iS2010_66a. 
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Fig. 17.— (upper-left) Frequency dependence of two sets of frequency differences, Au = 
j/obs _ j^mod^ between our fitted frequencies, and the corresponding Model S frequencies, 
j^mod, n = 0 ridge. The Au values computed using our Method 1 hts are shown as 

the black line, while the Au values computed using Method 2 are shown as the red line. 
The Method 1 frequency differences exhibit a pronounced discontinuity at the location of 
the switch-over in the htting range that was employed in Method 1. This location is marked 
by the vertical green dashed line in all three of the left-hand panels, (upper-right) Same 
as upper-left panel, but for the n = 1 ridge. In this panel, and in the other two right- 
hand panels, the vertical green dashed line marks the location of the transition from the 
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narrow fitting range to the wide fitting range in Method 1 for this ridge. In both top panels 
the dashed bine line indicates a freqnency difference of zero. Note that the vertical scales 
are different in the two top panels, (middle-left) Freqnency dependence of the linewidths 
compnted using Method 1 (black) and Method 2 (red) for the n = 0 ridge. The Method 1 
linewidths exhibit a sharp discontinuity precisely where the narrow htting range was replaced 
with the wide htting range. By contrast, the Method 2 linewidths do not exhibit any such 
discontinuity at the same frequency, (middle-right) Same as middle-left panel, but for the 
n = 1 ridge. The discontinuity in the linewidths occurs at the same frequency as in the 
upper-right panel, (lower-left) Frequency dependence of the logarithms of the amplitudes, 
or power densities, computed using Method 1 (black) and using Method 2 (red) for the n = 0 
ridge. The Method 1 amplitudes exhibit a discontinuity at the precise frequency where we 
switched from the narrow htting range to the wide htting range, (lower-right) Same as 
lower-left panel, but for the n = 1 ridge. The discontinuity in these amplitudes occurs at 
the same frequency as in the upper-right and middle-right panels. In all six panels, the hts 
were computed using the m-averaged spectral set iS2010_66a. 
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Fig. 18.— (upper-left) Radial order dependence of the normalized discontinuities in 
Method 1 (black diamonds) and Method 2 (red triangles) frequencies obtained from the m- 
averaged spectral set iS2010_66a. These normalized frequency discontinuities, / 2 )/F’ai/, 
were computed for the ridges n = 0 through n = 29 using the procedure that is described in 
Section Em The degree “/i” in all four panels is the degree of the highest-degree mode that 
we were able to £t for a given ridge using the narrow htting range in Method 1, while the 
degree “/ 2 ” represents the degree of the lowest-degree mode that we were able to fit for the 
same ridge using the wide fitting range in that method. In all four panels the two dashed 
green lines represent normalized frequency discontinuities of ±3, while the dashed blue line 
is plotted for normalized frequency discontinuities of zero, (upper-right) The same two sets 
of normalized frequency discontinuities that were shown in the upper-left panel are shown 
here as functions of the fractional inner turning-point radius, rt/i?©, of the li mode for each 
ridge. Here, Rq is the radius of the Sun. (lower-left) Same as upper-left panel, but for the 
normalized discontinuities in the linewidth w. (lower-right) Same as upper-left panel, but 
for the normalized discontinuities in the logarithm of the amplitude A. 
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Fig. 19.— (left) The frequency dependence of the logarithm of the Method 1 frequency 
uncertainties for the n = 1 ridge is shown as the black curve, while the frequency dependence 
of the logarithm of the corresponding Method 2 uncertainties is shown as the red curve. The 
vertical dashed line denotes the frequency where the wide fitting range replaced the narrow 
htting range in Method 1. (right) Same as the left panel, except that the two sets of frequency 
uncertainties were for the n = 2 ridge. Both sets of these frequency uncertainties resulted 
from fits to the m-averaged spectral set iS2010_66a. 
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Fig. 20.— Dispersion plane coverage of the set of frequencies obtained with Method 2 when 
applied to the m-averaged spectral set iS2010_66a. In the range 0 < / < 1000, 0 < n < 29, 
965 < f < 7000 /iHz we were able to successfully ht a total of 12,359 mode s . Th e black 
diamonds are for hts employing the asymmetric prohle of Nigam fc KosovichevI ( 1998 1. while 
the red diamonds are for those hts for which the symmetric Lorentzian prohle was used. The 
frequencies where the linewidths of most of the ridges exhibited either local or global maxima 
are shown as the green open circles that are connected by the green solid line. All of these 
frequencies can be seen to lie above the corresponding frequencies where the htting prohles 
were switched. The maxima in the linewidths are shown as a function of frequency in the 
upper-left panel of Figure |2T1 while the maximum in the n = 14 ridge is shown in the 
upper-right panel of Figure l22l 
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Fig. 21.— The frequency dependencies of the linewidths (upper-left), line asymmetries 
(lower-left), amplitudes, or power densities (upper-right), and powers (lower-right) of the 
12,359 hts that were made using Method 2 on the m-averaged spectral set iS2010_66a. The 
linewidths of many of the ridges exhibit either local or global maxima near frequencies around 
5800 to 6000 /iHz before decreasing slightly at even higher frequencies. For v < 4638 /iHz, 
only the asymmetric prohle was used and the vast majority of the line asymmetries were less 
than zero. For reasons described in Section 14.61 for 4638 < v < 4842 /iHz, the asymmetric 
prohle was used for some of the ridges, while the symmetric Lorentzian prohle was employed 
for others. For v > 4842 /iHz only the symmetric prohle was used, so B was set equal to 
zero for all of those cases. The dashed red line in the lower-left panel is for a line asymmetry 
B = 0. 
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Fig. 22.— (upper-left) Degree dependence of frequency, u, for the n = 14 ridge from fitting 
the m-averaged spectral set iS2010_66a using Method 2. The black diamonds represent the 
htted frequencies in the range 0 < / < 274. The red line is for a smooth £t to these htted 
frequencies, (lower-left) Degree dependence of the frequency slopes, Au/Al, for the same 
ridge. The values of Au/Al are plotted as the black diamonds at the abscissae / -|- 0.5, 
and are computed by subtracting two consecutive values of the htted frequency shown in 
the upper-left panel, i.e., Ai//A/|;+o .5 = i'{l -|- 1) — i^(/). The blue line is the derivative 
with respect to I of the smooth curve that has been htted to the frequencies and which is 
shown as the red line in the upper-left panel, (upper-right) Degree dependence of the modal 
linewidth, w. We note that the linewidth goes through a relative maximum at I ~ 167 
which corresponds to a frequency of 5686 /iHz. This is an example of the linewidth maxima 
that were illustrated in Figure [211 (lower-right) Degree dependence of the linewidth slopes, 
Aw/Al, plotted as the black diamonds at the abscissae I -|- 0.5. They were computed by 
subtracting two consecutive values of the htted linewidth, i.e., Atc/A/|i+o ,5 = w{l + l) —w{l). 
The blue line is the derivative of the smooth cnrve that has been htted to the linewidths 
and which is shown as the red line in the upper-right panel. 
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Fig. 23.— (left) Degree dependence of the differences between the observed frequencies, 
which were computed using Method 2 on the m-averaged spectral set iS2010_66a, and 
the theoretical model frequencies, in the sense of Sun minus model, for the n = 1 

ridge. The frequency differences, are shown as the black diamonds. The red 

line is for the Chebyshev polynomial, of degree w = 2A htted to the frequency 

differences. The 276 green triangles represent modes flagged as outliers. We note that these 
frequency differences are the same that were shown as the red curve in the upper-right panel 
of Figure [TTl (right) Same as left panel, but for the n = 11 ridge. For this ridge the degree of 
the htted Chebyshev polynomial, is -oj = 14, and 10 modes were hagged as outliers. 

In both panels the dashed blue line is for a frequency difference of zero. Note that the vertical 
scales are different on the two panels. 
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Fig. 24.— (upper-left) Degree dependence of the scaled frequency differences, as 

defined in equation (144)) . for the n = 1 ridge. The black diamonds are for the scaled frequency 
differences, while the red line is for the Chebyshev polynomial, of degree vr = 23 

that has been htted to the scaled frequency differences. In this panel all of the 435 outliers 
that were identihed in step 1 (276), and in steps 2 and 3 (an additional 159) are shown as 
the green triangles, (lower-left) Same as upper-left panel, but for the n = 11 ridge. For 
this ridge the degree of the htted Chebyshev polynomial, is vr = 14, and the 39 

outliers that were detected in step 1 (10), and in steps 2 and 3 (an additional 29) are shown 
as the green triangles, (upper-right) Degree dependence of the emulated scaled frequency 
differences, as dehned in equation ([59]). A comparison of the panels in the upper row 

indicates that the trendline of the emulated scaled frequency differences, as shown 

in the upper-right panel, is in excellent agreement with the trendline of the curve htted to 
the scaled frequency diherences, Sun^i, as shown in the upper-left panel. Therefore we can 
presume that the sharp variation of the scaled frequency diherences, dun^i, with degree at 
around / 60 is a real feature, (lower-right) Same as upper-right panel, but for the n = 11 

ridge. The march of the htted polynomial as shown in the lower-left panel is in excellent 
agreement with the march of the emulated scaled frequency diherences, shown in 

the lower-right panel. Hence, we can regard the sharp variations with degree of the scaled 












frequency differences, 5vn,h at about I ~ 1 and I ~ 7, respectively, as real features. 
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Fig. 25.— Location of the outlying cases in the l-u plane that have been detected by means 
of the outlier correction scheme as described in Section 17.11 The red diamonds are for the 
870 outliers for which both the frequencies and their uncertainties were adjusted and for the 
286 additional outliers for which only the frequencies were adjusted. The green diamonds 
are for the 1,129 outliers for which only the frequency uncertainties were adjusted and for 
the remaining 198 cases for which neither the frequencies nor their uncertainties needed to 
be adjusted. 
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Fi g. 26.— The relative squared soun d-speed deviations from the Standard Solar Model S 
of IChristensen-Dalsgaard et al.l (119961) as a function of fractional radius that we obtained 
by a structural inversion of the 6,366 frequencies and their uncertainties that resulted from 
the application of the outlier correction procedure to our original frequency table that has 
been computed using Method 2 on the m-averaged spectral set iS2010_66a, and covered the 
frequency range of 969 to 4500 /iHz. The horizontal bars represent the width (“spread”) 
of the localized averaging kernels, providing a characteristic of the spatial resolution; the 
vertical bars are the formal error estimates. 





















